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Nonperturbative Wilson coefficients of the operator product expansion (OPE) for the spin-0 glue- 
ball correlators are derived and analyzed. A systematic treatment of the direct instanton contri- 
butions is given, based on realistic instanton size distributions and renormalization at the operator 
scale. In the pseudoscalar channel, topological charge screening is identified as an additional source 
of (semi-) hard nonperturbative physics. The screening contributions are shown to be vital for con- 
sistency with the anomalous axial Ward identity, and previously encountered pathologies (positivity 
violations and the disappearance of the -+ glueball signal) are traced to their neglect. On the basis 
of the extended OPE, a comprehensive quantitative analysis of eight Borel-moment sum rules in 
both spin-0 glueball channels is then performed. The nonperturbative OPE coefficients turn out to 
be indispensable for consistent sum rules and for their reconciliation with the underlying low-energy 
theorems. The topological short-distance physics strongly affects the sum rule results and reveals 
a rather diverse pattern of glueball properties. New predictions for the spin-0 glueball masses and 
decay constants and an estimate of the scalar glueball width are given, and several implications for 
glueball structure and experimental glueball searches are discussed. 
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I. INTRODUCTION 



Among exotic hadrons, i.e. those which evade classification according to the constituent quark model, glueballs 
QHlll^l occupy an extreme position. In fact, they contain no valence quarks at all and are the only hadrons which would 
persist in a world without quarks. Even the exclusion of sea quarks from QCD, as in quenched lattice simulations, is 
not expected to alter their properties drastically [94[ • Hence glueball structure provides a unique source of information 
£SJ ' on nonperturbative gluon dynamics and may even shed light on the often elusive gluon component of the light classical 
hadrons (and potentially hybrids). A suggestive way to access this information is to investigate the role of known 
C*~) ■ coherent gluon fields, with instantons as the most prominent example, in glueball structure. 

The operator product expansion (OPE) 4] of glueball correlation functions provides an effective analytical frame- 
work for such investigations. In fact, it exhibits several exceptional features of the gluonium channels already at the 
Q_i qualitative level. The probably most instructive one is directly associated with the defining characteristic of the OPE, 
i.e. its factorization of the short-distance correlators into contributions from hard and soft field modes. Indeed, since 
the glueball's gluon content should be mostly nonperturbative, one might expect it to manifest itself primarily in the 
soft contributions, i.e. in the gluon condensates. Surprisingly, this is not the case: the condensate contributions are 
unusually weak H, IE an( i cannot fully reflect the nonperturbative nature of low- lying gluonia. This suggests that a 
. , major part of the nonperturbative physics is relatively hard and thus resides in the Wilson coefficients. The present 
paper contains a detailed analysis of such contributions and a comprehensive study of their impact on spin-0 glueball 
JH \ properties. 

Indications for the onset of nonperturbative physics in the scalar glueball correlator at unusually short distances, and 
the ensuing departure from asymptotic freedom, date back to the pioneering days of the QCD sum rule approach 0. 
The prototypical candidates for such hard, nonperturbative physics, direct instantons, describe tunneling processes 
which rearrange the QCD vacuum topology [lj} m localized space-time regions small enough to affect the x-dependence 
of the correlators over distances |z| <C Aqq D . Although this physics enters the Wilson coefficients, it was up to recently 
ignored in glueball (and other) sum rules, with the exception of an early explorative estimate in Ref . . 

The development of instanton "liquid" models 0] allowed several bulk properties of the instanton size distribution 
in the vacuum to be estimated, and the emerging scales were later supported by lattice simulations [Til IT^ . On the 
basis of these scales, the mentioned exploration of instanton contributions to a ++ glueball sum rule Q found them 
to be of considerable size and indicated their potential for improving the consistency with the underlying low-energy 
theorem. However, the involved approximations and especially the neglect of the crucial instanton-induced continuum 
contributions did not allow for reliable estimates of glueball properties. 

Only recently, the exact direct instanton contributions to the spin-0 glueball OPE (to leading order in K) were 
obtained and the instanton-induced continuum contributions derived [13j . On the basis of the resulting "instanton- 
improved OPE" (IOPE) the first quantitative QCD sum-rule analysis of direct instanton effects in scalar glueballs 
became possible. This analysis resolved long-standing consistency problems of the ++ sum rules with pure ly p er- 
turbative Wilson coefficients and led to new predictions for the scalar glueball mass and "decay constant" fs ^ n 
particular, the direct instanton contributions were found to increase the value of fs about threefold and supported 
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earlier indications for an exceptionally small scalar glueball size 

uAum The y also 

suggested a prominent role of 

instantons in the binding of the ++ glueball and generated scaling relations between its main properties and bulk 
features of the instanton size distribution 0] . A subsequent analysis of the related Gaussian sum rules , based 
on the same instanton contributions, confirmed some of these results and studied more detailed parametrizations of 
the phcnomenological correlator representation. 

Previous implementations of direct instanton contributions to hadron correlators, including those in the ++ glueball 
channel, relied on several standard approximations. However, there are reasons to suspect that these approximations 
may cause artefacts in the sum-rule results. We therefore provide a more thorough and systematic treatment of the 
direct-instanton sector below. Our subsequent comprehensive sum-rule analysis will indeed reveal a significant impact 
of these improvements on the predicted glueball properties. In addition, we will extend the study of nonperturbative 
Wilson coefficients to other glueball channels. While nonperturbative short-distance contributions are small in the 
tensor channel (mainly due to the absence of leading instantons corrections), they turn out to generate a rather 
complex pattern of new physics in the pseudoscalar channel. The analysis of this physics and the ensuing predictions 
for 0~ + glueball properties are further central objectives of our investigation and will provide additional insights into 
the role of (semi-) hard nonperturbative physics in glueball structure. 

The paper is structured as follows. In Sec. [H]we discuss pertinent general features of the glueball correlators and 
qualitative aspects of the instanton contributions. We also set up the Borel sum-rule framework. In the following 
Sec. IIIII we start our discussion of the IOPE by summarizing the known perturbative contributions to the Wilson 
coefficients. The next section contains a detailed derivation and analysis of the direct instanton contributions to 
both spin-0 glueball channels. Here we also implement for the first time realistic instanton size distributions and we 
explicitly renormalize the ensuing Wilson coefficients at the operator scale. The exceptional strength of the direct 
instanton contributions to the spin-0 glueball correlators provides an ideal testing ground for these improvements, 
which we expect to be useful in other hadron channels as well. (Previous work on direct-instanton contributions in 
the classical meson and baryon channels, and in particular the program starting with Ref. |l8j . treated all instantons 
as being of the same size.) 

In Sec. |?|we identify and implement new contributions to the IOPE coefficients in the h glueball correlator. 
As the instantons, these contributions are of nonperturbative origin and associated with the topology of the gluon 
fields. They arise from the screening of topological charge in the QCD vacuum and can be derived almost model- 
independently by requiring consistency with the axial anomaly |l7j . We present compelling evidence for the screening 
contributions to be an indispensable complement to the direct instanton contributions. 

In Sec. IVII we embark on a comprehensive quantitative analysis of the developed IOPE, including perturbative, 
direct-instanton and screening contributions. We start by discussing several characteristic features of its Borelmoments 
and their relevance for predicting glueball properties. Subsequently, we study the origin and role of the various 
subtraction constants (which appear in the lowest-moment sum rules) as well as their impact on the sum-rule analysis. 
We then match the IOPE Borel moments to their phenomenological counterparts, which include both perturbative 
and instanton contributions to the duality continuum and either one or two isolated, narrow resonances. The eight 
ensuing Borel sum rules are analyzed numerically, the predictions for ++ and h glueball properties are obtained 
and their significance for glueball structure and experimental glueball searches is discussed. Finally, in Sec. IVIII we 
summarize our main results and conclusions. 



II. GLUEBALL CORRELATION FUNCTIONS 



The QCD correlations functions in the main glueball channels arc 

IL G (x) = (0\TO G (x)O G (0)\0), (1) 

where the local, composite operators O g with G G {S,P,T} are the standard gluonic interpolating fields, i.e. those 
with the lowest mass dimension which carry the quantum numbers of the scalar (0 ++ ), pseudoscalar (0 h ) and tensor 
(2++) glueballs: 

O s (x) = as G% (x) G^ (x) , (2) 
Op( I )=« s G;(x)G^( S ), (3) 
O t (x) = 0£„ (x) . (4) 

Here 0" is the energy-momentum stress tensor of QCD and G is the dual of the (Minkowski) gluon field strength, 

G fiv — —E ^i/paG^ . (5) 
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Note that Os (Op) is proportional to the Yang-Mills action density (topological charge density 95]) of the gluon fields. 
The interpolators create vacuum disturbances with the appropriate quantum numbers whose propagation properties 
we will study below in the framework of a short-distance expansion. The factors of a s ensure renormalization group 
(RG) invariancc (at least to leading order in a s — g 2 / (47r) 96]). The corresponding Fourier transforms will be written 
as 

n G (Q 2 ) = i J d A xe l » x (0\TO G (x) G (0) |0) (6) 
with Q 2 = — q 2 . The above expressions (and the following ones, if not stated otherwise) refer to Minkowski space-time. 



A. Low-energy theorems 

The zero-momentum limits of both spin-0 glueball correlators are governed by low-energy theorems (LETs) which 
provide additional first-principle information and useful consistency checks for IOPE and sum rule analysis. In the 
scalar channel, the LET belongs to a class which was derived by Shifman, Vainshtein and Zakharov on the basis of 
renormalization group and scaling arguments [19| . These "dilatation" theorems apply to the zero- momentum limit of 
the correlators 

Ho (V) = * / d'xe^ (TO s (x) O (0)) npert , (7) 

where O is an arbitrary local color-singlet operator. The subscript refers to UV regularization by subtraction of the 
high-frequency (perturbative) contributions, i.e. 

1 t°° A 

U o (0) = - — Imn ( S )-Imn^ re '( S ) . (8) 

71" Jo S 1 J 

The associated low-energy theorems read 

n o ( ) = ^<o> (9) 

where do is the canonical dimension of O and bo = HN c /3 — 2Nf/3 is the lowest-order coefficient in the perturbative 
expansion of the QCD /3-function. 

This general class includes as a special case with O — Os the LET which governs the low-energy behavior of the 
scalar glueball correlator: 

Hs(Q 2 = 0) = '^(a s G 2 ). (10) 

The appearance of the gluon condensate and the sizeable factor in front render the numerical value of lis (Q = 0) 
both large and rather uncertain. The present range of values for the gluon condensate is about (a s G 2 ) ~ 0.25- 

0. 75 GeV 4 . With the value (a s G 2 ) = (0.07 ±0.01) GeV 4 used in 0, for example, one obtains an upper limit 
n s (0) ~ 11.17 (a s G 2 ) = 0.78 GeV 4 while n s (0) ~ 0.4 - 0.6 GeV 4 is probably more realistic. In the whole range of 
acceptable values, however, the low-energy theorem provides by far the largest soft contribution to the k = — 1 scalar 
glueball sum rule (see below). 

The zero-momentum limit of the pseudoscalar glueball correlator is likewise governed by a low-energy theorem, 
although of a rather different nature |9?J • It is usually stated for the zero-momentum limit of the correlator of the 
topological charge density 

Q(x) = -^G^G a ^ = ^0 P (x), (11) 

1. e. for the topological susceptibility 

Xt = i [ d 4 x (0\TQ (x) Q (0) |0) = -^U P (Q 2 = 0) . (12) 
J (8tt) 

In QCD with three light flavors and m Uj d *C m s , the low-energy theorem then reads |2lj ] 

Xt = U - ± ~ (qq) ■ (13) 

m u + md 
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This expression is based on the chiral Ward identities for the (anomalous) flavor-singlet and -octet axial currents and 
generalizes the classic large- N c result of Di Vecchia and Veneziano As a consequence, we have 

TIP (Q 2 = 0) = (8^) 2 {q q ) (14) 

which reduces to lip (0) = in the chiral limit. In Section IVll fand IYl CI in particular) we will analyze the consistency 
of IOPE and glueball sum-rule results with the LETs JTJJl and (jT4*|) . 



B. Qualitative aspects of instanton contributions 



Two properties of the spin-0 glueball correlators indicate on general grounds that they will receive relatively large 
direct instanton contributions: (i) the underlying spin-0 interpolators and couple exceptionally strongly to 
instantons and (ii) the leading-order (in K) instanton contributions are enhanced by inverse powers of the strong 
coupling (cf. Eq. I|65|)l. 

In addition, the instanton contributions show a strong dependence on spin and parity of the glueball correlator. 
This is reminiscent of the situation encountered in light meson and baryon correlators, where the remarkable topo- 
logical, chiral, flavor and spin-color structure of the quark zero-modes [22| in the instanton background generates a 
characteristic channel dependence pattern |23l l24j (which is difficult to reproduce even in more sophisticated quark 
models [2^]). Although the zero modes and quarks in general play a much smaller role in the glueball correlators, an 
equally (if not more) distinctive channel dependence exists among them, too. It is rooted in the (anti-) self-duality of 
the (anti-) instanton's field strength, 

GflV = iG^i; , (15) 

(in Euclidean space-time) and therefore of purely gluonic nature. (Anti-) self-dual gluon fields have color-electric and 
-magnetic fields of equal size, Ef — ±Bf . As a consequence of the Bianchi identity, they form a subset of all solutions 
of the Euclidean Yang-Mills equation. 

The impact of self-duality on the glueball correlators is twofold, as can be read off directly from the interpolating 
fields © - Q- First, during continuation to Minkowski space the Euclidean self-duality equation (|T5)l picks up a 
factor of 2, so that 

0^ (a?) = ±iC>{!^ (x) . (16) 

Therefore, direct instanton (or, more generally, self-dual) contributions to the pseudoscalar glueball correlator are 
equal in size and opposite in sign to those of the scalar correlator. This property provides a strong link between 
the IOPEs of both channels, with far-reaching consequences for our analysis. Moreover, self-duality implies that 
instanton-induced power corrections (condensates) to the IOPE of both ++ and h correlators are restricted to 
a few terms and cancel in their sum Double-counting of soft instanton physics in both Wilson coefficients and 
condensates of the IOPE is therefore excluded (with a few potential exceptions to be discussed below). Furthermore, 
since direct-instanton induced interactions turn out to be attractive in the scalar glueball channel, they must be 
repulsive in the pseudoscalar channel. 

Another obvious property of self-dual fields is that their energy, proportional to E 2 — B 2 , is zero (as expected for 
vacuum fields). Hence their energy-momentum tensor vanishes, i.e. 

T /J) (x) = 0. (17) 

This implies the absence of leading direct instanton contributions to the tensor glueball correlator. As a consequence, 
we have nothing to add to the standard sum-rule analyses (with solely perturbative Wilson coefficients) in the 
2 ++ channel 0, H||, except for a comment: the soft (i.e. condensate) contributions to the tensor correlator 
are conspicuously small, higher-dimensional power corrections vanish |27j and instanton-induced power corrections 
are absent |6J] . Nonperturbative contributions to the IOPE of this correlator should therefore primarily reside in the 
Wilson coefficients of higher-dimensional operators (starting with the gluon condensate) . However, direct instanton 
(and other hard, nonperturbative) contributions, although not forbidden by self-duality, are suppressed by minimally 
four factors of Aqcd/Q ~ 0.2 where Q ~ 1 GeV is the typical momentum scale in sum rule analyses. Nonperturbative 
contributions to the IOPE coefficients are therefore expected to be small, too, perhaps related to the relatively weak 
binding in the tensor channel observed on the lattice 28] . In fact, the hierarchy of interactions induced by self-dual 
gluon fields (i.e. attractive, absent and repulsive in the ++ , 2 ++ and h channels, respectively) agrees with the 
level ordering among the lowest-lying glueball states in the quenched lattice spectrum [TH |28| . 
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C. Spectral representation, Borel moments and sum rules 



In order to make contact with the glueball information in the IOPE, we match it to the spectral representation of 
the correlators, i.e. we construct sum rules. The spectral representation in the hadronic basis is conveniently written 
in the form of a dispersion relation, 

1 f °° , Imn G (-s) , , 

n G (Q 2 ) = -/ o ds^fy, (is) 

where Imn G (— s) is the imaginary part of the correlator at time-like momenta and where the necessary number of 
subtractions is implied but not written explicitly. Among the lowest-lying on-shell intermediate states which contribute 
to the spectral function p G (s) = ImII G (— s) /tt are the lightest glueball in the G channel (probably with admixtures 
of light quarkonium) and possibly light mesons with the same quantum numbers and sufficiently strong couplings to 
the gluonic interpolators. Beyond the isolated resonances, the multi-hadron continuum will set in at a threshold St- 
We therefore write the phenomenological representation of the spectral function as 

Im Yl { * h) {-s) = Im Ii% es) (s) + Im n£ ont) (-s) . (19) 

A simple and efficient description of most hadronic spectral functions, which is particularly suited for QCD sum rule 
applications (and beyond, e.g. for the parametrization and interpretation of lattice data |scLl3l| ). can be achieved with 
one or two narrow resonances and a duality continuum [29j. For the IOPE sum rules in the spin-0 glueball channels, 
we will limit ourselves to maximally two isolated resonances which we describe in zero-width approximation, i.e. as 
poles 

Imn% es) (s) = 7r/ G1 m G1 £ (a - m G1 ) + 7r/ G2 m G2 <5 (s - m 2 G2 ) (20) 

with (0|O G (0) \Gi (k)) = fGi m Gi- Invoking local parton-hadron duality, the hadronic continuum is replaced by its 
quark-gluon counterpart. (A different description of the resonance region, in terms of Goldstone-boson pairs, can be 
found in Ref. |32| . The resonance position is then roughly estimated by integrating the Goldstone-boson continuum 
up to the bound imposed by the low-energy theorem HlOj l.1 The duality continuum is obtained from the discontinuities 
of the IOPE by analytical continuation to timelike momenta, i.e. 

Im P4 c ° nt) {-a) = 6 (s - s ) Im ng OP£) (-a) . (21) 

The step function ensures that the continuum is restricted to the invariant-mass region where it is dual to higher-lying 
resonance and multi-hadron contributions. This duality region starts at the effective threshold sq. 

For accurate QCD sum rules it is generally insufficient to match phenomenological and IOPE representations of 
the correlators directly in momentum space. However, substantial improvement can be achieved by a Borel transform 



(in the convention of Ref. on both sides of the sum rules. The Borel transform improves the "convergence" of 

the IOPE, eliminates subtraction terms and, most importantly, implements an exponential continuum suppression 
which puts more emphasis on the glueball states. For the quantitative analysis of the spin-0 glueball channels one 
usually considers a family of sum rules, based on the Borel moments 

£ G , fc (r) =B[(-Q 2 ) fc n G (Q 2 )] (r), fc e {-1,0, 1,2} . (23) 

The selection of these four moments, which differ in their weighting of the intermediate-state spectrum, has practical 
reasons: higher powers k > 2 reduce the sum rule reliability (e.g. by shrinking the fiducial regions, see Section IVI Djl 
without yielding much further information, while k < — 1 would introduce additional subtraction constants which are 
not determined by the low-energy theorems. (The new subtraction term in the k = —2 sum rule of the ^ glueball 
channel is proportional to the derivative of the topological charge correlator at Q 2 = and therefore relevant for 
understanding the proton spin content |3^ (cf. Sec. IVI E|i The Borel moments satisfy the recurrence relation 

C G ,k+i (r) = -i-Ca,* (r) . (24) 

OT 
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By applying the operator in Eq. b'2'Mi to the dispersion relation (|18|) with the spectral function l|19(l . we obtain the 
phcnomcnological representation of the Borel moments as 

L G,k \ T , m Gi, jGi,so) - / G1 ro G1 e gi +/ G2 m G2 e 

-4-1^^(0) + - / ds^Imng^f-^e-^. (25) 

In order to determine the hadronic parameters m G i, fa and sq, one matches these moments to their IOPE counterparts 
C^,° k PE ^ (to be calculated below) in the fiducial r-region where both sides are expected to be reliable. This leads to 
the sum rules 

£gT E) M = C G h k (t; m Gl , f Gi , s ) . (26) 

After the standard renormalization group improvement, these sum rules are conveniently rewritten in terms of the 
continuum-subtracted Borel moments of the IOPE, 

AiOPE) 1 f'° , klm ^{IOPE) 



K Gik (T;so)^-6 k ,- 1 Tl^>(0) + - fc k Imn^^(- S )e- ST . (27) 

T Jo 

The subtraction constants IIq OPE ^ (0) are generated by the nonperturbative contributions to the Wilson coefficients 
(see below). In terms of the TZc,ki the glueball sum rules attain their final form 

K G<k (r; So ) - -4,-iIlg^ (0) + f G1 m G \ 2k e- m o^ + f^m^ e~ m ^ (28) 

which isolates the glueball contributions, the subtraction term for k = — 1 and possibly an additional resonance on 
the right hand side. 

The values of the subtraction constants n^^O) are determined by the low-energy theorems 1|1U|) and (|14|l . (The 
perturbative UV contributions to the IOPE are, as in Eq. © , removed by renormalization of the perturbative Wilson 
coefficients, cf. Section Till Al ) The large constant ITg (0) is known to dominate the scalar k = — 1 sum rule with 
purely perturbative Wilson coefficients [6j. This leads to a much slower decay with r than in the k > sum rules 
and thus to a much smaller glueball mass prediction (well below 1 GeV). The direct instanton contributions [l3| 
overcome this mutual inconsistency which had plagued the ++ glueball sum rules since their inception. We wil 
furtherl elaborate on the role of the subtraction constants in Section IYl CI 



III. IOPE 1: PERTURBATIVE WILSON COEFFICIENTS 



Our main tool for the QCD-based calculation of the spin-0 glueball correlators at short distances is the instanton- 
improved operator product expansion (IOPE). The general expression for the IOPE at large, spacelike momenta 
Q 2 = -q 2 > A QCD , 

n G (Q 2 )= J2 c { d G) {Q 2 -n)(d d ) , (29) 

d=0,4,... 11 

exhibits the characteristic factorization of contributions from hard and soft field modes at the operator renormalization 
scale fi: hard modes with momenta \k\ > fi contribute to the momentum-dependent Wilson coefficients Cd (Q 2 ;fi) 

while soft modes with \k\ < /i contribute to the vacuum expectation values ("condensates") of the operators Od of 
dimension d which are renormalized at /i. 

The perturbative contributions to the Wilson coefficients constitute the conventional OPE and will be discussed 
in the present section. In addition, the Wilson coefficients receive crucial nonperturbative contributions from direct 
instantons, i.e. those with sizes p < /i -1 [98j. Their evaluation is the subject of Section llVl Additional hard and 
nonperturbative contributions to the IOPE of the correlator, due to topological charge screening, are identified 
and evaluated in Section Ivl 



A. Perturbative OPE coefficients 



In this section we summarize what is known about the perturbative contributions to the Wilson coefficients 
C d S ' P ^ (Q 2 ) of the scalar and pseudoscalar glueball correlators. At present, perturbative contributions are avail- 
able up to maximally three loops, and for operators up to dimension d — 8. The accordingly truncated OPE has the 
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form 

n^Q 2 ) - + cf > (6<) + cP (o 6 ) + cP (6s) (so) 

where G = S, P denotes the (0 ++ , h ) glueball channel. 

Up to O (cq) (i.e. up to 3 loops - the two powers of a s from the interpolating fields are not counted here), the 
coefficient of the unit operator has the generic Q 2 dependence 



MG) _ 



ln(^ 



A 



A. 



(G) 



(31) 



were we have omitted irrelevant subtraction polynomials in Q 2 . Since the next-to-next-to-leading order (N 2 LO, i.e. 

O (ctg)) corrections are known only for the unit-operator coefficients Cq , strict compliance with the perturbative 
expansion would require their exclusion. Nevertheless, due to their exceptional size, the power-suppression of their 

~ (G) 

higher-dimensional counterparts C\ >4 and their special role in the duality continuum, they are usually taken into 
account. 

The radiative corrections to the gluon-condensate coefficient C4 are known up to O (a s ). The NLO corrections are 
of particular importance in 64 since the leading O (a®) term is momentum-independent and therefore contributes 
only to the lowest (k — —1) Borel- moment sum rule. The overall Q 2 dependence of C4 is 

-Q 2 ' 



C 



(G) 



;o 4 ) = Bjp 

The 3-gluon condensate term has the generic form 



Co + Cl In ( ^ 



(32) 



(33) 



while the 4-gluon condensates produce a term of the form 

C 8 (o 8 )=A)^. (34) 

Collecting all four contributions, we arrive at the general expression for the Q 2 -dependence of the OPE with 
perturbative Wilson-coefficients in both spin-0 glueball correlators, 



n^ c) (g 2 ) 



A - 


f Ax In 








OH 


- Ci In ( 











^ I ( ^ 



1 1 



+ In ( ^ 



(35) 



During renormalization group improvement in Sec. If 11 Bl the coefficients Ai — Di will acquire a logarithmic Q 2 - 
dependence due to the presence of the running coupling a s (and of the anomalous dimension 73^ in the case of the 
Ci). 

In view of potentially large quarkonium admixtures to physical glueball states we emphasize that the perturbative 
Wilson coefficients receive explicit quark loop contributions, both from radiative corrections and through the per- 
turbative /3-function during RG-improvement. Especially in the scalar channel ,99| the perturbative 7V/-dependence 
(Nf is the number of light quark flavors with m q < Aqcd) can be significant (see also Eq. (JTDJ). The soft quark 
contributions to the condensates and their TVy-dependence are very likely larger but more difficult to estimate. 



++ channel 



At present, the most accurate perturbative coefficients of the unit-operator, renormalized in the MS scheme and 
with threshold effects included, can be found in Refs. [l6Ll3§ |. For three light quark flavors (i.e. iV c = Nf — 3) they 
read 



A (S) _ 



(S) 



1 

bo 
4 " 



659 a s 
' ~36~T 

65.781- 



247.480 



4 S> = 



-10.1250 



(36) 
(37) 
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In addition to the rather recently calculated O (a 2 ) contributions, these expressions contain several corrections to the 
O (a s ) contributions which were implemented in older sum rules analyses. In Ref. [39| . the quark loop contributions 

were omitted (corresponding to Nf — 0) and contained a small error which was later corrected and 

implemented into 0. The O (a s ) contribution to A^f* was corrected once more in Ref. js^l- As a consequence, 
the numerical coefficient given above is about 25% larger than the one used in Ref. 0. 

The coefficients Bi of the lowest-dimensional nontrivial operator, i.e. the gluon condensate, receive important 
radiative corrections which were calculated in Ref. |39J for Nf = 0. Including the quark-loop contributions for Nf = 3 
they become 



B 



(S) 



4a s 



1 + 



175 a s 
36~ T 



<a s G 2 ) 



B 



(S) 



<a s G 2 ) 



(38) 



((a s G 2 ) = (c^G^G 12 ^)) As a consequence, the O (a s ) contribution to Bq is about 20% larger than the one used in 
Ref. 0. Note that the O (1) coefficient is Q 2 -independent. Hence, the leading power correction will enter all but the 
lowest Borel- moment (i.e. Ttg — l) sum rules solely via radiative corrections. 

The 3-gluon condensate term, again including O (a s ) radiative corrections, was first calculated for Nf = |39| and 
adapted to Nf = 3 in Ref. 16j, with the result 



cf^Sa^G 3 ), C[ S) =0. 



(39) 



(( 5 G 3 ) = (gfabcG^G^G^)) Both Refs. and use the Nf =0 value c[ S) = -58a 3 (gG 3 ) 39] instead. (Ref. 

otherwise sets Nf = 3.) Up to O (a s ), there is no explicit 7V/-dependence in Cq S \ 

The contributions from the highest-dimensional (d — 8) operators appear in the combination 

<a 2 G 4 ) s := 14((a s / afcc Gf ip G£ c ) 2 ) - ((a s /a fcc G^G^) 2 ) (40) 

of 4-gluon condensates and were calculated to leading order in the pioneering paper : 

= 87ra s <a 2 G 4 ) . (41) 

Estimates of the numerical value of (a 2 G 4 ) are customarily based on the vacuum factorization approximation 0, 



(42) 



The uncertainty in this estimate does barely affect the quantitative sum rule analysis since the impact of the coefficient 
(|4*TJl is almost negligible. 



2. + channel 



The lowest-order perturbative contributions to the unit-operator coefficient of the ^ correlator where first ob- 
tained in Re f. l5j . Subsequently, the calculation was extended to O (a s ) |36| and augmented by the 3- loop corrections 
[38| in Ref. |42l| (where also an error in the O (a s ) contribution was corrected). The result is 



-4, 



i(-r) 



1 + 20.75- 



305.95 



72.531- 



7T 

a. 



4 P, = 



-10.1250 



(43) 
(44) 



(Ref. [IBl contains a misprint of the sign of the ln 2 -term in its Eq. (2).) The leading-order contributions to agree, 
as expected, with those of the scalar glueball correlator. The NLO corrections of Ref. 38] affect Ao only weakly 
but increase Ai by almost a factor of four. The N 2 LO contribution increases Ag by about a factor of two. Finally, 
A2 ~ 10~ 3 for a /it ~ 0.1 and thus is practically negligible compared to Aq and A\. 

The gluon-condensate contribution was calculated to lowest order in Ref. [5j and up to O (a s ) in Ref. j3g, with 
the result 



B, 



47r-(a s G 2 ), 



B 



(P) 



9 



(45) 
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(The sign of the LO contribution to Bq is corrected and the NLO contribution is not determined by the renormalization 
group |3(|). 

The (<?G 3 ) term was calculated to lowest order in Ref. Q- The 2-loop contribution was first obtained for Nf = 
[36| and later for N f = 3 H2, with the result 

Cr = "8- 2 (v)V 3 ), Cf>=0. (46) 

(For Nf = one has C[ P) = 58tt 2 (a s /7r) 3 (gG 3 ) instead 0.) 
For the d = 8 contributions, which is associated with the combination 

(g 4 G 4 ) p = 2g i J abc f ade (Gl v G c a0 G d 'f u G e ' al3 + lOG^G^G^G^) (47) 

of 4-gluon condensates, one finds p| 

D { o P) = ~(9^) p , (48) 

which is (explicitly and to the given O (a s )) iV^-independent. (Ref. contains a typographical error in the definition 
of (g 4 G 4 ) p .) As in the scalar case, the condensate (g 4 G 4 ) p is conventionally factorized |5j according to 

({f abc G%G^) 2 ) ~ A (G® v G a ' t " / ) 2 , (49) 

(f bc G* Q G c ^) 2 ~ i (G^G°^) 2 , (50) 

so that 

(g 4 G 4 ) p = 30vr 2 (a s G 2 ) 2 (51) 

and thus finally 

= ^a a {a a G*f . (52) 

As in the scalar case, the uncertainty in Dq has practically no impact on the quantitative sum-rule analysis. 

We note in passing that the perturbative contributions to Cq increase in both spin-0 glueball correlators with the 
order of a s . The size ratios of the O (a°) : O (a*) : O (a 2 ) corrections to A Q P ^ are about 1 : 2 : 3, for example (at 
a s /n ~ 0.1). This type of behavior is often encountered in QCD perturbation series and might indicate the approach 
to the asymptotic regime. 



B. Borel moments, continuum subtraction and RG improvement 

The perturbative Wilson coefficients enter the glueball sum rules l|28[) through their contributions to the continuum- 
subtracted Borel moments <|27[) . In the present section we outline the calculation of these moments in the dispersive 
representation. The unit-operator coefficients contain the strongest UV divergences and therefore require the maximal 
number of subtractions, which turns out to be three (see below): 

n<r t} (<2 2 ) = n G (0) - Q 2 n' G (0) + \q 4 ^ (0) - ^ £° da ^ + ™ . (53) 

The further evaluation of the continuum contributions requires explicit expressions for the imaginary parts of the 
perturbative Wilson coefficients at timelike momenta. Those are obtained from the coefficients of the last section 
whose g-dependence resides in a combination of powers (g 2 ) and logarithms In (— g 2 //i 2 ) . The imaginary parts of 
the pure power corrections (i.e. the Bq, Cq and Do terms) are concentrated at s = and thus do not contribute to the 
duality continuum, while the logarithms acquire an imaginary part — n at time-like momenta. For s > we therefore 
find 



Iml4 pc) (-«) = -7T 



J \ \ [i* J J s 



A s 2 + 2A lS 2 h-i[ — + A 2 s 2 31n" — ) - ir 2 } + B 1 - -i 



(54) 
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The next step is to calculate the Borel moments of Hq C ^ (Q 2 ) according to Eq. I|23|) and to subtract the continuum 
contributions, which leads to 



i r°° 

£g (r; s ) = 4 pc » (r) - - / dss k lmH™ (s) <T" 



(55) 



After inserting the imaginary part, Eq. (JS2J, the perturbative continuum can be evaluated explicitly. The resulting, 
hnal expressions for Cc,k ( T ! so) are conventionally written in terms of the partial sums 



and the exponential integral 



Pk (x) 



E x (x) 



n=0 



dt- 



(56) 



(57) 



They can be found in Appendix 1X1 

Before implementing the continuum-subtracted Borel moments into the sum rules l|28|) , one still has to perform the 
standard renormalization group (RG) improvement. Using the leading-order perturbative /?- function, RG improve- 
ment amounts to the replacements [5J 



a -> a s (1/t) 



-An 



b In (A 2 r) ' 



(9G 3 ) 



as (1/r) 
" s (M 2 ) 



7/11 



(9G 3 ) 



(58) 
(59) 

(60) 



(Formally, the inclusion of NLO corrections to the /3-function would be more consistent for the evolution of the O (oq) 
contributions to the unit-operator coefficients. However, their impact is negligible in the kinematical region relevant 

for the sum rules.) Following the notation of Eq. i|27|) . the resulting Borel moments will be renamed Wq ^ (r, so). (It 
is common practice to perform the RG-improvement after the Borel transform, i.e. after the subtraction terms are 
eliminated. Interchanging this order leads to the same results, up to corrections of higher order in a s .) 

The final expressions for the continuum-subtracted Borel moments, as generated by the perturbative IOPE coeffi- 
cients, are (7 = 0.5772 is Euler's constant and the small contributions from the A 2 coefficients arc omitted): 



1. 



4.1(^o) 



[1 - Pl (sot)] + 2 -^P- [7 - 1 + £1 (sot) + e~™° + In (s t) Pl (s t)} 



T* T 

B (t) + B 1 (r) [7 + E x (s r)] - C Q (r) T + d (r) r 



7 -l 



TS 



Ex (s t) 



Do (t) 



(61) 



2. 



n Gfl ( T > S o) = [1 - P2 (SOT)] + — 

Si 



3 1 



— [1-po (s t)} + C - d [7 + #i (s t)] + D r 
r 



(62) 



3. 



<x ; (r, S0 ) 



6A 12Ai 
— 4 - i 1 - Pa (sot)\ + — — 

T T 

;P2 (sot) + ln (s t) p 3 (s t) 



~~F +El ( s ot) + Po (sqt) + -pi (s t) 
6 2 

-4 [1 - pi (sot)] + — [1 - po (sot)] - A) (63) 
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4. 



24A 



[1 - pi (s t)} 



48Ai 



7- 



25 
12 



+ E 1 ! (s r) + po (sqt) + \pi (s r) 



+ (s t) + -p 3 (s r) + In (s r) p 4 (sqt) 



2B 1 



[1 - P2 (SOT)] + 



°±[1- Pi (s t)]. (64) 



We conclude this section with a comment on the perturbative iVf-dependence and the choice Nj = 3 adopted 
above. At present, the identification of all experimental glueball candidates remains controversial 0, C| and lattice 
simulations including light quarks are still limited to relatively small lattices and poor statistics Hence, quenched 
lattice calculations [28} are currently the most reliable source for information on the glueball spectrum. One might 
therefore be tempted to argue that setting Nf = in the perturbative Wilson coefficients would allow for a more 
accurate comparison with the quenched lattice results. 

However, the values of all nonperturbative input parameters, i.e. the condensates, the subtraction constants and 
the bulk scales of the instanton size distribution, are deduced more or less directly from observables and therefore 
refer to the physical world with three active flavors. The iV^-dependence of all these quantities, furthermore, is not 
known well enough to allow for a reliable extrapolation to Nf = 0. The realistic case Nf = 3, besides being physically 
more relevant and allowing for quarkonium mixing effects, is therefore the only consistent choice for the perturbative 
Wilson coefficients. 



IV. IOPE 2: DIRECT INSTANTONS 



We now turn to our main objective, the analysis of nonperturbative contributions to t he W ilson coefficients. In 
the present section we study their prototype, direct instantons (with sizes p < p~ x ) [24[ |100| . which mediate fast 
tunneling rearrangements of the vacuum. Direct instanton contributions to the scalar glueball correlators and Borel 
sum rules were evaluated in Ref. and found to be instrumental in resolving mutual inconsistencies among the sum 
rules and with the underlying low-energy theorem. Moreover, they generate more stable and reliable sum rules, scaling 
relations between glueball and instanton properties and new predictions for the fundamental glueball properties. (See 
also subsequent work on the related Gaussian sum rules 

After a few introductory remarks on direct instantons and the semiclassical approximation, we will discuss pertinent 
aspects of the instanton size distribution n (p) in the QCD vacuum. In particular, we comment on the shortcomings 
of the "spike" approximation to n(p) which underlied all previous studies of direct instanton effects, and then im- 
prove upon it by constructing more realistic finite- width parametrizations which embody the essential features of the 
instanton size distribution as revealed by instanton liquid model and lattice simulations. 

Armed with realistic size distributions, we then embark on the calculation of direct instanton contributions to both 
spin-0 correlators, starting from x-space, and discuss their properties in some detail. In the course of this discussion 
we derive several new expressions for the instanton contributions. We also outline the derivation of the instanton- 
induced continuum contributions [13j and obtain the continuum-subtracted Borel moments for general instanton size 
distributions. Next we discuss, on the basis of some analytical results, generic effects of finite- width distributions and 
their impact on IOPE and sum rules. Finally, we renormalize the direct-instanton-induced Wilson coefficients at the 
operator scale p (which was ignored in previous studies of direct instanton effects). 

To start with, let us recall some basic features of the semiclassical approximation which underlies the evaluation 
of the direct instanton contributions. (For an introduction to direct instantons and more details see Ref. |24|.) To 
leading order in H, the (Euclidean) functional integrals representing the correlation functions are evaluated at the 
relevant action minima, i.e. the instanton and anti-instanton solutions closest to the space-time points linked by the 
correlator, with field strength (in nonsingular Lorentz gauge) |45| 

G$>° (*) = =^ (65) 

9s [{x-x a f+p^ 

(xq and p are position and size of the instanton, respectively, and r\ a ^v is the 't Hooft symbol H^)- The functional 
integrals then reduce to integrals over the collective coordinates of the saddle points, i.e. over xq^, p and the color 
orientation U. Quantum effects break the conformal symmetry of classical Yang-Mills theory and thereby generate 
a nontrivial measure n(p), i.e. the instanton size distribution, for the integral over p. The factor gj 1 in Eq. I|65|) 
enhances the direct instanton contributions to the correlators 10 by a factor a ,7 2 compared to the leading perturbative 
contributions. Since a s /7r ~ 0.1 at typical sum-rule momenta Q ~ 1 GeV, this enhancement is substantial and partially 
explains the important role of direct instantons in the spin-0 glueball channels. 
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The restriction of the semiclassical approximation to the saddle points associated with the nearest (anti-) instanton 
is justified by the hierarchy of IOPE and instanton scales. The distances \x\ ~ 0-2 — 0.3 fm -c ^-qcd 

accessible to the IOPE, in particular, are much smaller than the average instanton separation R ~ 1 fm (which is 
inferred from to the instanton density h I+ j ~ 1 fm -4 in the vacuum, see below). The relative diluteness of the 
instanton medium, characterized by the ratio p/R ~ 0.3, further reduces the impact of multi-instanton correlations. 
Multi-instanton corrections were calculated explicitly in the pseudoscalar meson channel (which strongly couples to 
instantons) in Ref. 0] and indeed found to be negligible. Nevertheless, we will discover a unique exception to this 
rule, generated by short-range correlations among topological charges (including those carried by instantons) in the 
QCD vacuum, in Sec. 



A. Instanton size distribution 



As alluded to above, the evaluation of direct instanton contributions to the IOPE coefficients requires the instanton 
(and anti-instanton) size distribution n I j{p) as an input. By definition, nj (p) specifies the average number Ni 
a i (p) dpd^xg of instantons with sizes between p and p + dp and positions between xq and xq + dxo in the vacuum, i.e. 

f A d " Nl tort 
dpd^xo 

Integrating over p, one obtains the instanton density 

711 = J Q dpni{p) = ^ ~ (67) 

where V4 is the space-time volume. Since the vacuum (in the thermodynamic limit) contains on average equal numbers 
of instantons and anti- instantons, the same relations hold for the anti-instanton size distribution nj{p) and density 
hj. We therefore omit the subscripts /, I and define 

n (p) = ni (p) = nj (p) , (68) 

which implies h = hi — hj. As a consequence, the total number of pseudoparticles (instantons and anti-instantons) 
in a vacuum volume V4 is 2hVi. 

In the IOPE context, the status of the size distribution n (p) is similar to that of the condensates: both are generated 
mostly by long-wavelength physics and characterize universal (i.e. hadron-channel independent) bulk properties of 
the vacuum. Hence the associated scales have to be imported from other sources (e.g. from phenomenology or the 
lattice). Due to their definition in terms of QCD amplitudes this is unambiguously possible, and due to their small 
number and universal character the IOPE's predictive power is only moderately reduced. 

The essential scales of n (p) are set by its two leading moments, i.e. the overall instanton density 



and the average instanton size 

1 

P = - 
n 



dpn (p) (69) 



dppn(p). (70) 



A large body of successful phenomenology in the instanton liquid model 01 has settled on the benchmark values 
p ~ 0.33 fm and h ~ 0.5 fm 4 which we will use throughout the paper. These scales imply a mean separation R ~ 1 
fm between instantons and are inside errors consistent with the still exploratory lattice determinations [Til Il2l | (see 
Ref. [43 for a critical assessment). 

The simplest parametrization which is able to embody these scales is the "spike" distribution n sp (p) = hS (p — p) |ij . 
Since the IOPE coefficients are not expected to be particularly sensitive to details of n (p) , the spike approximation 
has so far been exclusively relied upon in evaluating direct instanton contributions. Nevertheless, it is clearly an 
oversimplification and produces several artefacts, as we will see below. Moreover, a few features of n (p) beyond 
the moments n and p appear by now well enough settled to be implemented in more realistic parametrizations. 
Information on the behavior of n (p) at small p is available from perturbation theory in the instanton background 
[22I l5l| and implies the power-law behavior 

n(p)^V°- 5 (71) 
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(where bo = HN c /3—2Nf /3). Furthermore, lattice results and instanton vacuum model arguments [Tol . IT2I Esl Eil IHof 
indicate that the large-p tail of n (p) has a Gaussian falloff, 

n{p) p ^? Be- C ^tf. (72) 

Note, incidentally, that this strong suppression of large instantons (with sizes p ^> p) is in marked contrast to earlier 
dilute-gas assumptions [Hlj which proved inconsistent (i.e. IR- unstable) . 

An additional benefit of averaging the instanton contributions over a realistic size distribution is that it affords a 
simple and gauge-invariant way to implement the operator renormalization scale p, as we will discuss in Sec. IIVFI 
In the following subsections we are going to introduce two specific finite- width parametrizations. Both of them are 
uniquely determined by their behavior at small and large p and by their two leading moments n and p. This is a 
useful restriction since additional details of n (p) are neither reliably known nor should they have significant impact 
on the IOPE coefficients. 



1. Spike distribution 

Before discussing finite- width distributions, we mention a few obvious properties of the simplest approximation to 
the instanton size distribution, the zero-width or spike distribution 

n s P (p) = n6 (p-p), (73) 

which sets the sizes of all instantons equal to their average value p and has up to now been exclusively used in 
evaluating direct instanton contributions. Its two lowest moments determine this distribution uniquely. Moreover, it 
becomes exact in the large- N c limit of instanton liquid vacuum models. 

Most other features of the spike distribution, however, differ qualitatively from more realistic finite-width 
parametrizations: the median is equal to its mean, all higher (centralized) moments, including the variance 

(( P -pf)=0, (74) 

vanish identically, and the support both at small and large p is exactly zero. As we will see in Sections IIV El and IVII 
these properties induce several artefacts, including unphysical oscillations in the instanton-induced spectral functions 
at timclikc momenta. 



2. Exponential-tail distribution 

The minimal parametrization of a size distribution with power-law behavior at small p and an exponential decay 
at large p is 

"ex P (P) = ^f^J (£) n -P (- (f. + 1) I) ■ (75) 

All its additional properties and scales are determined by the two lowest moments (|69|l and (J7DJ. Although the 
exponential large-p tail is less favored by lattice and instanton liquid results than the Gaussian one to which we turn 
in the following section, it is nevertheless useful because it allows the p-integral in the instanton-induced spectral 
function to be done analytically. The resulting expressions provide a convenient benchmark for analyzing generic 
effects of finite- width distributions. 

For Nf = 3, instanton-background perturbation theory fixes the power of the small- p behavior at n = 4 (cf. Eq. 
GJ), i.e. 

n ^ = ¥3-p{f) 4 ^{- 5 j) i0rNf = 3 - (?6) 

The variance of the exponential distribution, 

((P - pf) = \J ncx P (p)p 2 dp - f = l -p\ (77) 
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implies a half-width of about 0.7 p, similar to lattice and instanton liquid results |l(J,ll2|. This is reassuring since the 
width of n cxp is not an independent parameter but determined by p. The maximum of the distribution is reached at 

Ppeak = "jT- (78) 

(i.e. Ppeak — 4/3 GeV -1 for the standard value p ~ 1/3 fm). A parametrization for n (p) with a similar exponential 
tail, 

N fa \ 
n(j>) = —expl---Pp), (79) 

has been used in a study of the asymptotic OPE behavior 35]. As our distribution JJSJ, it allows relevant p-integrals 
to be done analytically. However, it fails to reproduce the power-law behavior at small p and cannot be made to 
reproduce both n and p, thus introducing additional, insufficiently constrained parameters. 



3. Gaussian-tail distribution 



A size distribution with power behavior at small p and a Gaussian tail at large p can be obtained from the minimal 
ansatz 

n (p) = Np n exp (80) 

by requiring the two lowest moments to satisfy Eqs. I|69|) and (|70|l . This determines A and N uniquely, with the 
result 



n+2M»+ 1 



. 2 

p 



(81) 



Gaussian-tail distributions were found in the instanton liquid model (ILM) jlfj as well as in other approaches [Tol lisl 
|H3,|52. They are also supported by quenched lattice data [l^. In the ILM, in particular, one finds 

n ILM (p) - Cp b °- 5 exp £T) (82) 



2 f 

(bo = llA c /3 — 27V//3) which approaches the spike approximation at large TV, 



At N c — Nf = 3 it becomes 



lim niLM [p) oc 5 (p - p) . (83) 



to be compared with our parametrization l|81|) for n = 4, 



. . 2 18 n /p\ 4 / 2 6 p 2 . 

(P) = -SaZs = 7 ex P (iV> = 3) . (85) 



3 6 7T 3 /3 \p J \ 3 2 7T p 2 

The width of the ILM distribution is about 10% larger. Equation (|85|l is the most realistic parametrization among 
those which we will consider. 

The distribution (|85|) has the variance 

3 2 5tt A 1 



15 



which corresponds to a half-width of about p/2, somewhat smaller than the half-width of the exponential-tail distri- 
bution l|7t)|) . The position of the peak lies at 

3 fn 

Ppeak = -J-p~0.Mp. (87) 

Since n g is more symmetric than rt oxp , its p pea k is somewhat larger. 

Due to the stronger decay at large p, n g (p) contains a larger number of intermediate-size instantons with p ~ p. 
Since only instantons with p < p^ 1 ~ 2p contribute to the IOPE coefficients (p, is the operator renormalization scale), 
the Gaussian-tail distribution - cut off at p^ 1 - will typically result in larger instanton contributions to the Borel 
moments than the exponential one. 



B. Direct instanton contributions in a>space 



In this section we derive, to leading order in the semiclassical expansion, the direct instanton contributions to the 

~ (s P) 

IOPE coefficients Cq ' of the spin-0 glueball correlators. We also comment on contributions of higher orders in 
H and on coefficients of higher-dimensional operators. The calculations, as well as the discussion of the singularity 
structure and of the pointlike-instanton limit, are best performed in Euclidean space-time, the instanton's natural 
habitat. 

As outlined above, the semiclassical contributions from instanton-induced saddle points to the functional integral of 
the glueball correlators are obtained by evaluating their Wick expansion in the instanton background, e.g. by means 
of the gluon background-field propagator of Ref . [53j . Adding up the contributions from the nearest instanton and 
anti-instanton one finds, to leading order in h, 



II ' (x 2 ) =Y.J d P n (P) J <^ x * ( T °s (*) O s (0)) J+/ (88) 



A^ 1 ) f..:2 

i J 
2 9 3 2 



dpn(p) / d 4 x Q - \ — . (89) 

[(x-xo) + p 2 f[xl + p 2 ] 4 

We recall that it suffices to perform all calculations in the scalar channel since the pseudoscalar correlator receives, 
due to self-duality, identical instanton contributions (in Euclidean space-time) . A convenient analytical expression for 
the above xo-integral in terms of a hypergeometric function |54| . 

o 12 7T 2 / 9 

J \{x-Xq) - 

is derived in Appendix 151 



d X °U ^2 , our 2 2H = Z? 2i M 4 ' 6 '9'-Z^ ' ^ 90 ) 

[(x-x ) + p 2 ] 4 [xl + p 2 ) 4 42 V 2 W 



The singularity structure of the instanton contributions at fixed instanton size p can then be read off from 



by recalling the well-known analyticity properties of the hypergeometric functions |54j : Eq. (|91|l has no singularities 
inside a circle of radius 2p around x = 0, but it has cuts emerging from the two branch points x = ±2ip. (These 
singularities have their origin in the (partially gauge-dependent) denominators of the instanton's field strength |J?jSJ|.) 
Contributions from these branch points dominate the Fourier transform of Eq. I|91|l at large momentum transfer (see 
next section). The final result for the instanton contributions is obtained by integrating Eq. (|91|l over all instanton 
sizes with the weight n (p) : 



nf f) (* 2 ) - J dpn ip) !!('+'-) (x 2 ;p) = ^jdp n -f >Fx ( 4 , 6, £, - £ 



(92) 



Note that the p-integration over physically reasonable instanton size distributions, with the power-law behavior 
()71JI at small p and Nf < 3, does not produce singularities of Eq. I|92|l at x 2 = 0. This implies, in particular, that the 
x 2 — > limit is finite (see below). From 2F1 (a, b, c, 0) = 1 one obtains 
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The ^-dependence of Eq. (|91|l shows that the instanton-induced correlations are maximal at x 2 — and strongly decay 
(oc \x\~ at fixed p) for \x\ ^> 2p, i.e. when the arguments of the correlator cease to lie both within the bulk of the 
instanton's localized action density. Note, furthermore, that the x-dependence of the fixed-size instanton contribution 
reduces to a delta function in the small- p limit, 



limn 



(/+/) , 2 v 2 8 3 1 



•<4 - 



(94) 



This pointlike-instanton approximation will be useful below. 

Technically, the result (|92|1 provides the leading nonperturbative contributions to the IOPE coefficient Cq S ^ (x) of 
the unit-operator in the scalar glueball channel. (Strictly speaking, soft-mode contributions, in particular those due 

to large instantons with p > pT x , should still be excluded, cf. Sec. IIV Fl ) Of course, the coefficients of higher- 

~ (S) 

dimensional operators receive direct instanton contributions as well. The coefficient C\ of the lowest-dimensional 
condensate (aG 2 ), for example, gets a correction from the process in which one of the gluon fields emanating from the 
interpolator J5J) is soft and turns into a gluon condensate while the other one is hard and propagates in the instanton 
background. The general form of this contribution is 



n s!<aG 2 ) ( x2 ) i a sG 2 ) I dpn(p) I d 4 x 



(x-x ) + p 2 



(95) 



Relative to the unit-operator term, this correction is suppressed by four powers of e = \x\ Aqcd ^0.2. More generally, 
power corrections associated with d-dimensional operators are parametrically suppressed by a factor e d and a lso b y 
the relatively large glueball mass scale (cf. Sec. IVlfl . Therefore, we neglect such contributions in the following |101| . 

The IOPE coefficients also receive instanton corrections from higher orders in the semiclassical expansion. To O (h), 
these corrections arise from Gaussian quantum fluctuations around the instanton fields |lfj| and have the form 



n 



(/+/) 

G,0(h) 



(x 2 ) 



G#'« (x) [DPDVD va (x,y)] ab G^ b (y) 



(96) 



where is the covariant derivative in the adjoint representation and D v p is the gluon propagator in the instanton 
background. Radiative corrections of this sort are suppressed by at least two powers of a s . Moreover, the average 
instanton is small on the QCD scale, p~Aqcd <C 1, and the instanton action Si (p) ~ 10ft consequently large compared 
to the action of the quantum corrections. We therefore do not expect such corrections to have significant impact on 
our analysis and do not consider them further. 

Finally, it is instructive to compare our result i|92[) with direct instanton contributions to pseudoscalar meson 
correlators |46l l55l |56j . The latter arise mostly from quarks propagating in the zero-modes states of the covariant 
Dirac operator in the instanton background 22]. (The non-zero modes can normally be approximated by plane 
waves.) In the meson sector, the distinct topological, chiral, flavor and spin-color properties of these zero modes 
are responsible for most features of the direct instanton contributions, including their strong channel dependence. 
Clearly, the instanton-induced physics in the glueball correlators is quite different in this respect. (It does not require, 
incidentally, to account for ambient soft vacuum fields in the gluon propagators which, in contrast, rather strongly 
affect the quark zero-mode propagation.) 



C. Borel moments 



The next step in our program is to calculate the Borel moments associated with the direct instanton contri- 
butions (|92|l . This is best done by first transforming to (Euclidean) momentum space, 

n 5 +/) (Q 2 ) = J d 4 xe lQa; 4 /+/) (x 2 ) (97) 

or, more explicitly, 

= 2\ 2 [ dpn(p)(Qp) 4 K 2 (Qp), (99) 
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where J\ (z) and K2 (z) are Bessel and McDonald functions 54] . In the last equation we have made use of the integral 
153 



f°° / 9 r 2 \ 7 

J drr 2 J 1 {Qr) 2 F 1 [^-,--A = _QV #f (Qp) . 



(100) 



(An alternative calculation of Eq. I|99|) , which agrees with the expression first obtained in [|| , can be found in Appendix 

El) 

It is instructive to examine the limits of (Q 2 )- Using for simplicity the spike distribution (|73[l and the 

asymptotic behavior of the McDonald functions [54| , one finds 



4 /+J) (Q 2 ) ^ 2Vn {Qpf e~ 2Q P. (101) 

The exponential decay of the integrand at large Q is expected since II^. + ' (x 2 ; p) is analytic in the circle \x\ < 2p 
around x 2 = 0. Its scale is set by the singularities nearest to the real axis, i.e. the branch points at x = ±2ip. 
Physically, this can be understood as the phase-space suppression encountered when distributing the hard momentum 
Q over multiple soft modes of the coherent instanton field with size p. While the fixed-size instanton contribution 
1|91JI contains no asymptotic power corrections (due to the absence of singularities at x 2 =0), integration over the 
instanton size does produce inverse powers of Q. However, they start with Q~ 5 (the power is determined by the 
large- p behavior H71|) of n (p)) and therefore do not interfere with the power corrections of the truncated IOPE (cf. 
Eq. ES|). 

In the opposite limit, i.e. for Q 2 —> 0, the instanton contributions turn into 

n s /+7) (Q 2 = 0) = 2V J dpn (p) . (102) 

This term plays a key role, as an instanton-induced subtraction constant, in the lowest (i.e. k = —1) Borel moments 
and in the associated spin-0 sum rules (cf. Sec. IVI C(l . The obligatory removal of soft contributions from instanton- 
induced Wilson coefficients (cf. Sees. HVFI and lVI Cfl de-emphasizes contributions from multi-instantons and other soft 
vacuum fields (not taken into account in Eq. (|92|l ') and should thereby improve the nearest-instanton approximation. 
Although the latter is strictly valid only for \x\ <C R, it often seems to work over larger distances (as long as cluster 
decomposit ion does not become an issue [H|)), due to the strong localization and small packing fraction of the 
instantons |102| . Nevertheless, one should not expect Eq. I|102H to be a complete representation of the correlator at 
Q = 0, as can be seen, for example, from the one-instanton approximation to the gluon condensate, 

(aG 2 ) I+I = 16tt J dpn (p) , (103) 

which, when multiplied by 327r/&o, yields only about half of Eq. (|102fl and therefore does not saturate the low-energy 
theorem (|10H . The inconsistency with the low-energy theorem (I14fl in the pseudoscalar channel, incidentally, is far 
more dramatic and can be overcome only by additional nonperturbative physics, cf. Sees. Ivl and lYl CI 
According to the general definition l|23|l . the Borel moments of the direct instanton contributions are 



-Q 2 f4 I+I) (Q 2 



(104) 



(for k G { — 1,0,1,2}). In order to calculate these moments explicitly, it is convenient to rewrite Eq. (|99|l in terms 
of an integral representation for the McDonald function before applying the Borel operator (|22|l . In Appendix El we 
outline this calculation for k = — 1, which results in ^| 



4£? (r) = -2V/ dpn{p)ee-t 



(:i+0/v'„(0+(2 + e + |)A"i(0 



(105) 



Above we have introduced the dimensionless variable £ = p 2 /2r. From the lowest, "generating" Borel moment all 
higher moments are obtained by differentiation with respect to — r, 



4'tll W = "J^f } (r) (k > -1), (106) 



and explicitly read 
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For r — > 0, at fixed p, all instanton-induced Borel moments vanish exponentially. This is a direct reflection of the 
exponential large-Q suppression in Eq. i|101|) and renders these moments practically negligible for r < 0.2 GeV~ 2 (at 
p = p). More specifically, if £ = p 2 /2r ^> 1 and if n (p) is replaced by the spike approximation {73J, Eqs. I|105|) - 
((109(1 reduce to 
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(110) 



In the case of k = 0, for instance, Eq. I|110fl is reliable when r < 0.4 GeV~ 2 but it underestimates the maximum at 
t ~ 0.6 GeV -2 and has the wrong decay power at larger r (about twice the correct one, see below). The estimate 
of the instanton contributions to the k = — 1 sum rule in Ref . [9( were based on the approximation 1110(1 . Due to 
the wrong decay behavior and the missing continuum contributions, however, it is unsuitable for reliable glueball 
mass and coupling estimates. In fact, after subtracting the associated continuum contributions (as in Eq. 1)55(1. see 
next section), the exponential suppression l(110|) will turn into an enhancement and generate increasing instanton 
contributions down to r = 0. 

At intermediate r (r > 0.2 GeV~ 2 ), (r) rises quite steeply towards its maximum. The associated, fast 

variations of its slope generate large oscillations in the higher Borel moments (due to the r-derivative in Eq. 1(104(1 ') 
which would increasingly impede reliable sum rule fits. Again, this problem is overcome by subtracting the crucial 
instanton-induced continuum contributions, which transforms the oscillations into a monotonic decay. In the opposite 



reaches the negative 



limit, i.e. for r p /2, the k > Borel moments approach zero from above while Cg _ 1 

subtraction constant — + (0). The general decay behavior of the instanton contributions for large r is an inverse 
power law. This continues to hold for fairly general n (p) (including ((76(1 and 185(0 as long asr> p 2 /2 where p c is 
roughly the maximal size at which n (p) has appreciable support. More explicitly, one has 
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with positive coefficients afc. As a consequence, the instanton-induced Borel moments decay as fast as their leading 
perturbative counterparts for r — > oo, and faster than all power corrections (cf. Section IlII B(l . (Of course, this limit 
is not physical since the short-distance expansion breaks down for r — > oo.) 



D. Imaginary part and instanton-induced continuum 

The discontinuities of the IOPE at time- like momenta generate the duality continuum ((21|) . In the spin-0 glueball 
channels, the instanton-induced continuum contributions were first obtained in Ref. |l3;j and shown to have crucial 
impact on consistency and predictions of the ++ -glueball sum rules. In the following, we outline their calculation 
and analyze some of their pertinent properties. Besides providing new insight into the structure of the instanton 
contributions, this analysis will prove useful in assessing the effects of realistic instanton size distributions in subsequent 
sections. Finally, we assemble the instanton-induced, continuum-subtracted Borel moments which enter the spin-0 
glueball sum rules. 

The imaginary part of the instanton contributions is obtained by analytically continuing Eq. ((99(1 in the complex 
s = q 2 = —Q 2 + ie plane. The behavior of the McDonald functions under analytical continuation can be expressed as 


I =^e- ivv l 2 Hi l) (ze-™/ 2 ) for f < argz < ir V ' 
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where Hu is the Hankel function of the first kind, 

HP (z) = J v (z) + iY v (z) 
From the cut structure of the Hankel functions [54[ one then finds 
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which leads to [l3| 
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Several qualitative and quantitative features can be readily deduced from this expression. First of all, for s 
(i.e. s <C p~ 2 ) it reduces to 
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which has the leading s 2 -dependence of the free gluon loop. This implies, for example, that the large-r behavior of 
the corresponding, continuum-subtracted Borel moments (see below) equals that of the leading perturbative Wilson 
coefficient. Since the continuum subtraction does not affect the falloff at large r, it also confirms the decay behavior 
HI lip established in the last section. 

A more important property to notice is that the analytical continuation has turned the exponentially small instanton 
contributions (|99|l at large Q 2 and fixed p into strong oscillations with increasing amplitude at timelike s. This becomes 
more explicit at large s, 



Imn( /+/ ) (-a) ^? 2 4 7r 3 S 3 / 2 np 3 cos {2y/ip) 



(117) 



where we have used the spike approximation (|73|l . Such a behavior is familiar from the analytical continuation of 
semiclassical tunneling amplitudes in quan tum mechanics and can lead to a selective enhancement of parts of the 
amplitude (when crossing Stokes lines) |58j . Potential problems of this sort (including the strong oscillations) will be 
tamed by finite- width instanton size distributions (cf. Sec. II V F|) . In any case, the overall growth of the instanton- 
induced imaginary part oc s 3 / 2 at large s (and fixed p) is weaker than that of the free gluon loop (cf. Eq. (|54[) ) and 
the dispersive integral exists without subtractions, i.e. 
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The analytical evaluation of the s-intcgral indeed reproduces Eq. (|99|l . 

As discussed above, the small-Q 2 limit of the correlator 199|) lies outside the range of validity of both the nearest- 
instanton approximation and the IOPE. Usually, this does not cause concern because the sum rules are analyzed at 
momenta Q 2 > 1 GeV 2 where both approximations should work. However, the lowest Borel moment (with k = —1) 
by design contains the Q 2 = limit of the correlator as a subtraction constant. This term is not eliminated by the 
Borel transform since the latter is applied to n (Q 2 ) /Q 2 '- 
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which of course also follows directly from Eq. (|118|l with B 
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(cf. Eq. i|118|) at Q 2 — 0) and in the opposite limit 

Hm (r) = (Q 2 = 0) (123) 

(in agreement with Eq. (|lllfl ). 

Note that the zero-momentum limit H102(l of the correlator can be recovered from the dispersive representation 
(|TT£)) as 



where we have made use of the integral |57j 



dpn(p), (124) 



daaJiU/ap)Y i {y/ap) = j. (125) 

In order to avoid subtraction constants of this type, moments with negative k are usually not considered in sum-rule 
applications. In the spin-0 glueball channels, however, the k = — 1 moments are particularly useful because the 
subtraction terms on the phenomenological side (cf. Eq. (J2HJ) are determined by the low-energy theorems of Sec. 
Ill Al As shown above, the implementation of this additional first-principle information comes at the price of dealing 
with the IOPE- induced subtraction constants, too. Their interpretation, treatment and role in the k = — 1 sum rules 
will be the subject of Sec. IVI CI 

The dispersive representation of the higher Borel moments follows from Eq. I|121(l by taking the appropriate number 
of r-derivatives, according to the recursion ({106(1 . This results in 

4s t ^ ( r ) = -4,-in(; J+/) (Q 2 = 0) + - (°° dss k Iml4 /+/) (-s) e- ST (126) 

Jo 
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dpn (p) - 2 4 tt 3 J dpn (p) p A J dss k+2 J 2 (y/ip) Y 2 (y/ip) e- ST . (127) 
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Note that the UV-convergence of the dispersion integrals is ensured by the Laplace kernel and that Eq. I|127|) is 
often easier to handle than the individual integrated expressions (|105|) - (I109|l . Moreover, the subtraction of the 
instanton-induced continuum (cf. Eq. I|27|l ). 

(t^o) = 4T ] W - - [°° dss k lmlli 1+l) (s)e-»T, (128) 

provides now immediately a compact expression for the instanton contributions to all ++ glueball sum rules: 

^fc +7) ( T ; so) = -2V<J fc ,_ 1 J dpn (p) - 2 4 tt 3 J dpn (p) p 4 J ° dss k+2 J 2 {yfsp) Y 2 (^p) e - ST . (129) 

(These contributions are not affected by perturbative (one-loop) RG-improvement since the anomalous dimension of 
the interpolating field Os (x) vanishes to this order.) 

Depending on the choice for n (p) , the size of the instanton-induced Borel moments (|129|) is either significantly 
larger or comparable to the size of their perturbative counterparts (|61l) - (I64|l . The instanton continuum turns out to 
be indispensable for the sum rules to match at intermediate and small r [13| since the perturbative continuum cannot 
smoothly extend the exponentially vanishing instanton contributions l|127|) towards r — > (cf. Eq. ip.iup ). 



E. Analytical results for specific instanton size distributions 

In order to learn more about the impact of finite- width distributions n (p) on the spectral function (|115f) and the 
derived Borel moments, it will prove useful to carry out the integration over the instanton size analytically. As 
mentioned in Sec. IIV Al this becomes possible when specializing to the exponential-tail distribution (|76|l . First, 
however, we are going to establish the point of reference for later comparison by briefly reviewing the results for the 
spike distribution J73J). Inserting it into the instanton-induced spectral function l|115|l results in 

Im n L + P ? (~ s ) = ~^ 4 np 4 s 2 J 2 (V~sp) Y 2 (V~sp) , (130) 
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which was used in [13( and subsequent work. The plot of this expression in Fig. 1 (a) shows, besides the expected 
rather violent oscillations with increasing amplitude (note the scale), a strong violation of positivity for s > 4 GeV 2 . 
(We use the standard instanton scales n = 0.75 x 10 -3 GeV~ 4 and p = 1.69 GeV -1 throughout the paper.) At small 
s, the spike-induced imaginary part starts out quadratically, 

In.n£?h)-2Vn(p 2 S ) 2 . (131) 

At large s the rise becomes slightly weaker and gets modulated by a harmonic oscillation, cf. Eq. 1)117(1 . The direct 
instanton contributions to the continuum-subtracted Borel moments simplify to 



n st ^ ( T ) = -2Vn4.-i - 2\ 3 np 4 f ° dss k+2 J 2 (yfip) Y 2 (V«p) e- ST . (132) 
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Now we derive the analogous results on the basis of the more realistic finite-width distribution 1(76(1 with exponential 
large-p tail. Although a Gaussian falloff is favored by instanton-liquid and lattice simulations (cf. Sec. IIV A(l and 
modifies the detailed behavior of the instanton-induced imaginary part at s < p~ 2 , the resulting expressions will 
be a useful benchmark for assessing qualitative effects of realistic n (p) (and for evaluating them numerically). The 
exponential-tail distribution 1(76(1 specializes the imaginary part (|115|) to 

Imn^J (-a) = - 2 -^^s 2 I"" dpp*J 2 (VSp) Y 2 (VS P ) e-^'K (133) 
6 P Jo 

A somewhat lengthy but essentially straightforward calculation shows that the p-integral can be expressed as a 
combination of three hypergeometric functions, 



2 ' 2 ' ' Asp 2 J V 2 2 Asp 



(134) 



This spectral function is plotted in Fig. 1 (b). Comparison with its counterpart (|130|l in Fig. 1 (a) reveals that the 
finite-width distribution has turned the problematic, spike-induced oscillations into a monotonic falloff at large s. 
More specifically, owing to 2 F\ (a, 6, c, 0) = 1, the first and leading term in 1(134(1 decays as 

T tt( /+/ )/ % 3 3 5 7 7 2 11tt 4 n _ s/2 

ImII s,cxp (-*) — W 5 i5 s ' ■ ( 135 ) 

z p 

The damping of the large-s oscillations is a generic effect of finite-width distributions. Nevertheless, the imaginary 
part 1(1 34|) still changes sign once, due to the interference among the individual terms. The associated positivity 
violations are much milder than those in Eq. 1(130(1 . however, and their impact on the Borel moments in the relevant 
r-region is strongly reduced. Moreover, the imaginary part of the full IOPE is now nonnegative for all s, as it should 
be. We therefore expect that extended size distributions will improve the consistency of the Borel sum rules. 
At small s, Eq. ((134(1 becomes 

t rr( /+/ ) l \ 273 ' 771-3 - f-2 \ 2 

ImII S,c X p (~ s ) — ► p n (P s ) ■ ( 136 ) 

It is instructive to compare the above behavior to the s — > limit ((131(1 of the spike distribution: the quadratic 
s-dependence is maintained but Eq. I(136D is about three times larger. This enhancement can be traced to the large-p 
tail of rt xp (p) by expressing the small-s behavior of the instanton-induced spectral function for general n (/?) as 

Imn^ +/) (-s) s< ^ 2Vn (p 4 ) s 2 , (137) 

where ( p 4 ) is the fourth moment of the instanton size distribution. Comparison with Eq. 1(131(1 or direct evaluation 
via Eq. ((76|) shows that 

which relates the enhancement factor to the increased weight of instantons with p > p. We therefore anticipate this 
factor to be reduced when the large-p contributions are excluded from the IOPE coefficients in Sec. IIVFI 
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In view of the rather dramatic enhancement and modulation of the instanton-induced spectral function at large 
s by the spike distribution, one might wonder how the latter can produce useful approximations to the continuum- 
subtracted Borel moments (I129|l at all. The reason is twofold: first, these moments are essentially Laplace transforms 
of the spectral function, 

Ksk^ ( r ; s o) = -2V4,-i f dpn (p)-~ f ° dss k Imll^^ (-s) e" ST , (139) 



which implies that contributions from the s ^> t _1 region are exponentially suppressed. In addition, the large-s 
contributions of the spike- induced spectral function partially cancel each other due to the modulating oscillations. 
The impact of the large-s behavior on the Borel moments is therefore strongly reduced (except at very small r which 
are irrelevant for the sum rule analysis, see below). 

In order to explore the effects of finite- width distributions on the r- and so-dependence of the Borel moments more 
quantitatively, we again resort to the exponential-tail distribution which yields 
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Obviously, all Borel moments reach finite limits for r — > 0. Towards r — * oo we find from Eq. I|137f) that the 
r-dependence of the Borel moments vanishes as 

^s,t l] (r, -so) ™ -2 7 7r 2 4,-i | (p) + 2 3 tt 2 (fc + 2)!^J [1 - p fc+2 (, r)] , (141) 

where we have used the definition (|56|l for the partial sums pk (x). Of course this result holds for all reasonable 
instanton size distributions: specific choices for n (p) do not affect the r-dependence but just the overall magnitude 
of the moments at r ^> p 2 . 

The impact of finite- width distributions on the so-dependence of the Borel moments (which generates the ± -channel 
dependence of the direct instanton contributions) will play an important role in our subsequent discussion. From Eq. 
(I140[) we find that in the sum-rule relevant (r, so)-region (specifically for so > 2 — 4 GeV 2 and r > 0.2 GeV -1 ) the 
magnitude of the instanton contributions decreases more strongly with increasing so if the exponential-tail distribution 
is employed. The direct instanton contributions to the pseudoscalar glueball sum rule, where sq is typically a factor 
of two larger than in the scalar channel, therefore become smaller than those to the scalar sum rule. This reduces the 
instanton-induced repulsion in the ^ channel without significantly weakening the important attraction in the ++ 
channel. 

Furthermore, we observe that for sufficiently large sq the instanton-induced Borel moments l|140|) turn negative at 
small r. (For the k = moment this happens if r < r c ~ 0.3 GeV -2 . For k > 0, t c increases due to the stronger 
s fe -weight of the spectral function at larger s.) This positivity violation (which grows with sq and therefore has more 
impact on the pseudoscalar sum rule) would reduce the fiducial r-domain of the associated Borel sum rules (cf. Sec. 
(|VID|) ) but does not a priori prohibit a sum rule analysis. Nevertheless, it exposes a shortcoming of Eq. <|140|) which 
can be traced to the contributions from large instantons. In fact, we will argue in the next section on general grounds 
that such contributions have to be excluded from the IOPE coefficients. (The discussion of the far more serious 
instanton-induced positivity violations at large t in the ^ channel will be the subject of Sec. IVI ) 

Finally, we note that the above results remain qualitatively unchanged when n cxp (p) is replaced by the Gaussian- 
tail distribution (|85|l . Even the quantitative results are generally very similar, as expected from the arguments of Sec. 
II V Al and confirmed by numerical integration of the corresponding p-integrals. 



F. Implementation of the operator renormalization scale 



In the exact version of the (I)OPE, contributions to the Wilson coefficients originate exclusively from hard field 
modes with momenta larger than the operator renormalization scale p, |J,|59J. Potential double-counting of soft modes 
and IR renormalons [6(| are thereby excluded from the outset and the //-independence of (RG-invariant) short-distance 
amplitudes becomes manifest. However, the gauge-invariant implementation of IR cutoffs in the perturbative Wilson 
coeff icients is technically complex and commonly neglected in QCD sum-rule applications (with a few exceptions 
103]). This "pragmatic" approximation works well since in QCD soft perturbative amplitudes are generally much 
smaller than their nonperturbative counterparts (which generate the condensates). 
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However, it is a priori unclear whether the pragmatic neglect of IR cutoffs works for nonperturbative Wilson 
coefficients as well, as tacitly assumed in all previous work on direct instantons. In the following we will argue that 
realistic, finite-width instanton size distributions permit an explicit and straightforward implementation of such IR 
cutoffs for the direct instanton contributions. This allows us to improve upon the "pragmatic" treatment of the 
nonperturbative sector and to assess its conceptual basis and range of validity. (The exceptional strength of the 
direct instanton contributions makes the spin-0 glueball correlators a particularly suitable testing ground for this 
renomalization. ) 

Our implementation of the renormalization scale p is based on the observation that an instanton with fixed size 
p mainly contributes field modes with momenta |fc| < p^ 1 to the correlator. The contributions of large instantons 
with p > p~ x and the fluctuations around them are therefore almost exclusively soft and should be excluded from 
the Wilson coefficients. This (partial, see below) cutoff procedure is gauge invariant simply because the instanton 
size is. Since the renormalization scale associated with the conventional condensate values (cf. Sec. IVI All is not very 
accurately determined, it would not make sense to treat /i _1 as a sharp cutoff. Instead, we implement it smoothly by 
replacing the full instanton size distribution n (p) with 

^(p) = fj (p- 1 -p)n(p) (142) 

where the soft step function 9p can be chosen, e.g., in the form of a Fermi distribution 
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The "diffuseness" parameter (3 sets the scale for the width of the transition region. With (3 <C p the cutoff practically 
ceases to exist, i.e. all instanton sizes are about evenly affected, whereas for (3 ^> p a sharp cutoff is reached. For 
practical calculations we find values (3 ~ 5/5" 1 ~ 3 GeV to be an effective compromise between these extremes. (For 
(3 > 5 GeV strong oscillations of the imaginary part set in (even at s > 20 GeV 2 ) while for (3 < 2 GeV the cutoff is 
already rather ineffective.) Alternatively, we will use 

Op (p- 1 - p) = 1 - i arctan [(3 (p - p,' 1 )] , (144) 

which has a softer transition region and a bigger large-/? tail for equal values of (3, to estimate the sensitivity of the 
sum-rule results to the details of the cutoff procedure. 

The physical interpretation of the above /i-implementation is very transparent: the modified distribution 1142|) just 
exclu des non-direct instanton contributions to the IOPE coefficients since those are already included in the condensates 
|104| . Nevert heless , this does not exclude all \k\ < p modes since even arbitrarily small instantons contain some soft 
contributions 105]. A complete renormalization would naively amount to replacing 

(*)->/ ^e-^&JHk) (145) 

where G/P is the Fourier transform of the instanton field. However, such a procedure is unacceptable since it violates 
gauge invariance. In fact, devising a gauge-invariant cutoff is a complex task since gauge transformations generally 
change the momentum content of the fields. This problem is well-known in the application of Wilson's renormalization 
group to gauge fields and has not yet been solved satisfactorily. Nevertheless, there is reason to believe that our above, 
simplified renormalization of the direct instanton contributions is sufficiently complete since the overall momentum 
transfer Q acts as an additional IR cutoff by suppressing field modes with momenta k <C Q. 

An interesting question in this context concerns the direct-instanton induced /(-dependence. The /(-dependence of 
both the perturbative coefficients and the c onde nsates (which arise from the anomalous-dimensions) is logarithmic 
(cf. Appendix and therefore rather weak jl06l |. Since mutual cancellations among all contributions are required to 
render the correlators /(-independent, one expects the /t-dependence of the nonperturbative coefficients to be similarly 
weak. This turns out to be indeed the case, mainly because large-/? contributions are already strongly suppressed by 
the external momentum scale Q > 1 GeV in Eq. I|99(l and by the exponential tail of realistic instanton size distributions 
(cf. Eq. G2J0- Hence the p^ 1 cutoff does only reduce intermediate-size instanton contributions significantly, implying 
a relatively weak /(-dependence for p < pr 1 (which holds for p = 0.5 GeV and pT 1 ~ 0.6 GeV adopted below). 

One could contemplate, incidentally, to identify p, = Q in Eq. 11421) by analogy with the RG improvement of 
the perturbative IOPE coefficients. However, avoiding or summing large logarithms is not an issue here, the impact 
on the sum rules would be limited by the moderate /(-dependence, and the additional s-dependence in the spectral 
integrals of Eq. (|129|l would obscure the comparison with spike distribution results. We will therefore keep p fixed in 
the quantitative analysis below. 
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The implementation of p has several significant effects. One of them is the reduction of the instanton-induced 
imaginary part for s <C p 2 (where the excluded soft-mode contributions had mainly contributed), and a second one 
is the shift of its leading peak towards larger s. Both features can be understood by noting that 



Imlli /+r) (s) s<< ^r 2 2 3 7rW 



for s <C fM , where we have defined the moments of the renormalized size distribution n (p) as 
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The reduced magnitude of the imaginary part (|146|) is obviously due to Op (/Lt 1 — p) < 1 for all p, which implies 
\p k ) < \P/- The shift of the peak is caused by the stronger reduction of the higher moments: due to the asymmetry 

of n (p), impressed by its limiting behavior (cf. Eqs. H71fl an d (|72|) '). <^p fe ) jp k increases with increasing k. Since this 
enhancement originates from the large-p tail of the distribution, it is strongly reduced by renormalization and the 
ratio (p k )^ I (p k ) therefore decreases with k. 

Another significant renormalization effect is the reduction of the number of "active" instantons in the Wilson 
coefficients, i.e. 



dpn (p) 



dpn^ (p) = n d ir = (n 
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with £ < 1, where fidir correctly accounts for the density of direct instantons. This effect is missed when using the 
spike distribution since it is normalized to the total instanton density, i.e. large instantons with p > p are not excluded 
but simply counted as instantons of average size. (This suggest to use fidir instead of n to improve direct instanton 
calculations on the basis of the spike distribution.) The above reasoning indicates, incidentally, that direct instanton 
calculations involving the spike distribution strictly make sense only for p < [/, . Although satisfied by the (not very 
accurately determined) standard values p~ x ~ 0.4 fm and p ~ 0.33 fm, this condition would be violated by larger 
choices for p and/or p. 

Summarizing the main results of this section, we find the exclusion of non-direct instantons by an explicit large-p 
cutoff to be mandatory when implementing finite-width distributions. However, the "pragmatic" neglect of p together 
with the use of the spike distribution (with p < p -1 ) can be useful for an approximate estimate of the Borel moments 
since large instantons (with p 3> p) are rare in the QCD vacuum and since their contributions are bounded by 
Q -1 acting as an additional cutoff. This justifies the tacit assumption underlying earlier implementations of direct 
instantons. Nevertheless, the renormalization generates several new and important effects, including an effective 
reduction of the instanton density, which will significantly improve both sum rule consistency and results (see Sec. 

ED. 



IOPE 3: TOPOLOGICAL CHARGE SCREENING 



The direct instanton contributions to ++ and <~ glueball correlators differ only in their sign (cf. Sees. Ill Bl 
and I1V|) . As a consequence, the instanton-induced attraction in the scalar channel turns into a strong repulsion in 
the pseudoscalar channel. Therefore it is hardly surprising that the signal for the pseudoscalar glueball disappears 
and that even the general positivity requirement is violated |42j when these repulsive contributions are added to the 
perturbative Wilson coefficients. We demonstrate this by plotting the k = Borel moment in the h channel, 
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obtained from both perturbative and instanton-induced Wilson coefficients (on the basis of the spike distribution 
in Fig. 2 (a). The sum-rule relevant so-range starts at so > 6 — 8 GeV 2 , where the Borel moment indeed 
becomes negative for r > 0.8 GeV -2 , i.e. in the middle of the fiducial domain (cf. Sec. IVI|) . (The negative areas 
of the Borel moment at smaller sq and r are an artefact of the spike distribution which will be removed by realistic 
instanton size distributions, cf. Sec. I VII ) Moreover, the r-slope becomes positive for r > 1 GeV -2 , and both features 
make a match to the decaying exponential of a resonance (cf. Eq. (|28|l ') impossible. Contrary to lattice results |28l |43| . 
the corresponding sum rule would therefore predict the absence of low-lying pseudoscalar glueballs. In addition, we 
will find below that the low-energy theorem i|14fl in the h channel would be badly violated. 
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In search for the origin of these problems one should keep in mind that the (up to the sign) identical expressions for 
the direct-instanton induced Borcl moments l|129|) . involving the same approximations, play a highly beneficial role 
in the ++ sum rules and enhance their consistency I n fact, this might not be as paradoxical as it first appears: 
although the moments l|129|) have the same magnitude in both channels, the direct instanton contributions to the sum 
rules do not. The difference lies in the value of the duality threshold sq which is typically more than twice as large in 
the h than in the ++ channel. (This is mainly a reflection of the larger h glueball mass.) Their rather strong 
so-dependence thus equips the direct instanton contributions with a certain channel dependence. 

As demonstrated above, however, this channel dependence is too weak to resolve the mentioned deficiencies. Yet it 
is conceivable that the underlying approximations, i.e. the spike distribution and the neglected renormalization, may 
have underestimated the so-dependence. Partly for this reason we have improved upon these standard approximations, 
by implementing realistic instanton size distributions and the renormalization scale. It turns out that their combined 
impact indeed mitigates the problems in the ^ glueball channel, but it is too small to resolve any of them. 

Thus we face for the first time the situation that direct instanton contributions seem to worsen the consistency 
of a set of QCD sum rules and to create serio us n ew problems. Since our treatment of the direct instanton sector 
leaves not much room for further improvement |107| , one is bound to conclude that essential additional physics is still 
missing in the IOPE coefficients of the glueball correlator. Our main guide in the search for this physics, which 
will be subject of the next section, is the fact that it apparently did not show up in other hadron channels: previously 
considered sum rules were found to be consistent and stable with perturbative and direct-instanton induced Wilson 
coefficients only. 



A. Topological charge correlations 

As argued above, the failure of direct instanton contributions to generate acceptable pseudoscalar glueball sum rules 
strongly suggests that equally important - and therefore very likely nonperturbative - contributions to the Wilson 
coefficients are still amiss. Since these contributions should affect almost exclusively the h glueball correlator, one 
is led to search for a unique property which singles this channel out among other hadron correlators. The result is 
quite obvious: the interpolating field Op is proportional to the topological charge density Q, 

Op (x) = 8irQ(x) , (150) 

which implies that the h glueball correlator is proportional to the topological charge correlator |I08| . Hence it 
distinguishes in the strongest possible way between instantons and anti-instantons and is highly sensitive to instanton- 
anti- instanton correlations. (In contrast, Os (x) is proportional to the gluonic action density and therefore treats 
instantons and anti-instantons on the same footing.) 

The correlations between topological charges are generated by light quark exchange or, equivalently at low energy, 
by the exchange of flavor-singlet pseudoscalar mesons. The form of these interactions is dictated by the axial U (I ) 
anomaly 66] 

^2 9i7^759i = 2 ^2 miqiij 5 qi + 2N f Q. (151) 

i—u : d,s i=u 1 d,s 

Compliance with the ensuing behavior under Ua (1) transformations requires the lightest pseudoscalar flavor-singlet 
meson tjq to couple to the topological charge density as [l7| 

AC = -i-y^Q (x) % (as) (152) 

where the coupling 7 ?(0 is often written as 7 I)0 = ^2Nf / f Vg . The form of this interaction and most of its qualitative 
consequences do not depend on the specific origin of Q (x) . Instanton models and lattice simulations indicate that the 
topological charge density in the QCD vacuum is rather strongly localized and mostly carried by instantons [Icl H3| j 
however, and we will adopt this picture in some of the developments below. Owing to the coupling i|152|) . ?7o exchange 
generates strong correlations between topological charges which add to the direct-instanton induced ones (of range p) 
evaluated in Sec. IIVI The quark chiralities entering the anomaly equation (|f 51 \ require the ?7o-mediated interactions 
to be attractive (repulsive) between topological charges of opposite (equal) sign. 

Between isolated topological charges these interactions would be of long range since the 770 would be a light (quasi- 
) Goldstone boson. In the topological-charge ensemble of the QCD vacuum, however, the 770-exchange forces are 
screened by the formation of Debye-Huckel clouds |68j | in which positive (negative) topological charges surround 
themselves with negative (positive) ones. This collec tive mechanism, in turn, renders the ?7n-induced interactions 
short-ranged by generating a screening mass for the 770 |l09j | (which solves the U (1) problem (69j). Topological charge 
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screening has been discussed in Refs. [73, E3 an d investigated in ILM [T^ as well as unquenched lattice [T^, 
simulations. Indirect evidence for strongly localized screening clouds has also been found in the diffcrcnce-to-sum 
ratio of the isovector scalar and pseudoscalar ILM correlators where "unquenching" strongly suppresses the 
instanton-induced spin- flip probability for a propagating quark at short distances. 

In order to understand how topological charge screening affects the IOPE, one has to compare the factorization 
scale p with the characteristic scale of the screening mechanism. Since the physical rj' mass (i.e. the 770 mass after 
r/o — r]s mixing is taken into account, see below) is about twice as large as p, topological charge screening takes place 
over distances of the order of the screening length 

A scr ~ — 0.2 fm < p- 1 (153) 

which are small on the IOPE scale. Hence screening contributes mainly to the Wilson coefficients. This also fol- 
lows directly from the observation that Q-screening strongly influences the x-dependence of the topological charge 
correlator at |x| < A scr <§; ^-qcd^ which resides exclusively in the Wilson coefficients. And since Q-screening is a 
nonperturbative and collective mechanism originating from the light-quark sector, it obviously cannot be described 
either by perturbative or by direct instanton contributions. Thus we have reached our objective to uncover nonpertur- 
bative contributions to the Wilson coefficients which go beyond the direct instanton approximation and affect almost 
exclusively the h glueball correlator. 

Under the premise that the topological vacuum charg e den sity is mostly associated with instantons, the screening 
corrections can be identified with multi- instanton effects |110| . In fact, topological charge screening has been observed 
in instanton vacuum model calculations and found to affect mostly the h glueball and rf channels [ll], ^| . Even if 
light-quark induced interactions and therefore screening effects are neglected, phenomenological and lattice results in 
the majority of all calculated hadron correlators are well reproduced. However, in the 7/ and the related h glueball 
channel (and to a lesser degree in the scalar-isovector do meson channel) this approximation clearly fails and the rf 
correlator turns negative. "Unquenching" of the light- fermion sector resolves these problems and generates the 77' |7q 
and h glueball [15| resonances (although neglected interference terms between classical and one-loop effects did not 
allow to determine h glueball properties in Ref. [l5|). 



B. Screening contributions 



The screening contributions to the topological charge correlator (and therefore to the + glueball correlator) can 
be calculated from the pertinent low-energy approximation to the (Euclidean) QCD generating functional 



Z[0]= J D{A^q,q,C}exp(^-S Q cD + t I 



d 4 x9 (x) Q (x) 



(154) 



in the presence of a classical source 9 (x) . The correlation function of the topological charge density (and the related 
one for the h glueball) is then obtained by functional differentiation, 



(Q(x)Q(O)} =-Z[0}- 1 



s 2 z[e] 



66 (x) 89 (0) 



(155) 



8=0 



We now outline the derivation of the appropriate low-energy approximation to Z [9] . (More details can be found 
111 Refs. mini-) To b egin with, one separates the domain of the gluonic integration variable into the underlying 
multi-instanton background and the non-topological remainder, over which one integrates first. (Of course, in practice 
this could be done at best approximately, on the lattice, but the part of the result which is relevant in our context will 
be obtained indirectly from the Ua (1) anomaly l|151|) and effective field theory arguments.) It remains to sum over 
all classical multi-instanton configurations and to integrate over their collective coordinates (U, p, Xo) i . Integrating 
over the color orientations Ui and sizes pi turns the i-th (anti-) instanton into a colorless lump of topological charge 
density centered at xqj. Equation Ijl54(l then becomes the partition function of a medium of localized topological 
charges (which we approximate as pointlikc), interacting via a multi-local "potential" V: 
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d 4 x ,i exp -V (x Qi i...x 0;N++ N-) +i^2Qi9{x s) , 



(156) 



where n is the average density of either positive or negative topological charges. In order to find an approximate 
expression for V, we recall that the lowest-lying excitations of the QCD vacuum are the pseudoscalar Goldstone 
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bosons. Since (in the chiral limit) the flavor-singlet 770 gets its mass from the topological gluon fields over which we 
did not integrate yet, we expect it at this stage to be part of a degenerate U (3) flavor nonet. At low energies, the 
interactions V between the topological charges will then be dominantly mediated by 770 exchange, 

V (x 0tl ...x ,N + +N-) - ^2 v Vo (x ,i,X 0t j) , (157) 

i<j 

with the coupling to Q (x) determined by the axial anomaly (cf. Eq. (|152[l 'l. Since the 770 contains a superposition of 
both quark chiralities with opposite sign, the anomaly implies that equal (opposite) topological charges repel (attract) 
each other. We therefore have 



, d q 7 

v no {x ,i,x j) = qiQj I 7—4 2 — — 2 — exp (iq (x ,i - xqj)) (158) 



2 

Vo 

(2tt)* q z + m' NG 

where qi — ±1 is the topological charge located at xo,i and m 2 NG is the part of the 770 mass due to finite current 
quark masses. After absorbing an infinite factor into the renormalization of n and performing the remaining sums 
and integrals in Eq. (|156|) . one ends up with the low-energy approximation 



and the effective action [77J 

Seff [110,6} = I (fx 



N / D770 exp (-S e// [770,0]) (159) 



(160) 



Expanding the action i|160|) up to O (77^) one finds that the interaction between the topological charges in the vacuum 
has generated the expected screening contribution 

m 2 scr = 2ntf (161) 

to the 770 mass, m^ = 77i^ G + mj. cr , which does not vanish in the chiral limit and solves the U (1) problem. (The tree 
approximation in Eq. I|158|l is reliable as long as fv)t <C 1 and n/m% G ^> 1.) With the help of Eq. (|155|) we finally 
obtain the low-energy approximation 

(Q (x) Q (0)) = ^ 27^ (x) - (2n 7 „o) 2 (Vo (x) Vo (0)) (162) 

(8tt) 

to the topological charge correlator. Although we have referred to instantons as the source of the topological charge 
density Q (x) for definiteness, the above arguments would hold with inessential modifications (e.g. in the set of 
collective coordinates) for other localized, topological-charge carrying gluon fields as well. 

The first term in Eq. (|162[> is just the correlation due to the topological charge "cloud" of a single instanton in the 
pointlike limit (which we have adopted above for simplicity). This can be seen explicitly from the direct instanton 
contribution 192fl by implementing the pointlike-instanton limit (|94fl with the help of the spike distribution for p = 0, 
i.e. n p i (p) = n5 (p). The result 



inf (8tt) 2 



dpn p i (p) U y p > [x 2 ; p) = 2n<5 4 (x) (163) 



is indeed identical to the first term in Eq. i|162[l . The second term in Eq. (|162fl is the expected Debye screening 
correction. 



C. Singlet-octet mixing, Borel moments and subtraction term 

The screening contributions to the topological charge and h glueball correlators in Eq. (|162fl do not yet contain 
the 7/0 — 778 mixing corrections which arise from finite light-quark masses. These corrections are rather substantial, 
however, and should be implemented into the IOPE. This will also allow us to use the experimental 77' and 77 meson 
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masses and mixing angle, instead of extrapolations to the chiral limit. In addition, it is more consistent with our use 
of the standard con dens ate values below which were obtained from phenomenology and thus correspond to realistic 
light quark masses Hi3. The quark -mass dependence of the perturbative Wilson coefficients, on the other hand, 
originates from (dominantly strange) quark-loop corrections and should be negligible. 

The screening contributions to Eq. (|163|l can be adapted to physical quark masses by expressing them in terms of 
the 77 and rf mass eigenstates. These are related to the singlet (770) and octet (77s) flavor eigenstates by 



|v7o) \ = /cosy - sin tp\ ( \rf) 
W) J V s[n( P COS( P J V Iv) 



(164) 



where tp is the 770 — 77s mixing angle. The 770 correlator therefore becomes 

(770 (x) Tjo (0)) = cos 2 <p (rf (x) 77' (0)) + sin 2 <p (77 (x) 77 (0)) , (165) 
which allows us to rewrite the screening contributions to the topological charge correlator (|163fl in the form 

(Q (x) Q (0)) scr = - (2n 7t)n ) 2 [cos 2 <pD [m v ,,x) + sin 2 tpD (m v ,x)] , (166) 
where we have used the Euclidean tree-level propagator 

(77 (*) 77 (0)) = D (m n ,x) = -^K x (m n x) (167) 

for 77 and the analogous one for 77'. Diagonalization of the physical Goldstone boson mass matrix [toL relates the 
screening mass (|161H to the physical masses of the pseudoscalar mesons as 

mlcr = 2 nJno = m %' + m % ~ 2m K- ( 168 ) 

(Note that Eq. I|168|) reduces to m 2 cr = m 2 = m 2 / in the chiral limit.) After continuation to Minkowski space-time, 
our final expression for the screening contribution to the IOPE of the pseudoscalar glueball correlator becomes 

Tl ( p cr) (x) = F*,D (m v >,x) + F^D (m v , x) (169) 

where {F n i,F v ) — IGirnj^g (cos ip, sin (p) 

All mass and coupling parameters in Eq. (|169|l will be fixed at their experimental values. For the 77 and rf m asses 
and the mixing angle we use m, = 547.30 ± 0.12 MeV, m,- = 957.78 ± 0.14 MeV and ip = 22.0° ± 1.2° 78]. The 
axial anomaly renders even the Q-770 coupling and the overall coupling 2ti7^ experimentally accessible, with the result 
(2nj Vo ) 2 = 9. 732 x 10~ 4 GeV 6 [7^|. Alternatively, this coupling could be estimated from the standard value n = 0.5 
fm~ 4 = 7. 53 x 10 -4 GeV 4 of the instanton density and the experimental pseudoscalar meson masses according to 

(2n lvQ ) 2 = 2n (m 2 n , + - 2m 2 K ) , (170) 

which yields (with m K o = 497.67± 0.31 MeV Q) the value {2n^ m) ) 2 = 1. 086 x 10~ 3 GeV 6 . It is reassuring that both 
estimates are perfectly consistent. We will use the rounded value (2^) = 10~ 3 GeV 6 in our quantitative analysis 
below, which implies F 2 = 0.543 GeV 6 and F 2 = 0.0886 GeV 6 . 
From the Fourier transform of the screening contribution, 



4 scr) (Q) = 7W^ + ?5r^' < I7I > 



we finally obtain the Borel moments 

n^ ) (r) = S k ,- 1 f^ + Sl +F^e- m >+F^ m fe-^. ,172, 



The r-independent term in Eq. (|172|l is the screening-induced subtraction constant —IV < 'p Cr ^ (0) which will play a 
pivotal rule in the 72.p,_i sum rule analysis below. Already at this stage, however, its necessity can be seen from yet 
another inconsistency which one would encounter by restricting the nonperturbative Wilson coefficients exclusively 
to direct instanton contributions. Indeed, the latter generate a large subtraction constant (cf. Eq. I|102|) ') 

n£ +/) (0) = -2Vn (173) 
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(where we have adopted the spike distribution) which cannot be matched on the phenomenological side of the sum rule 
since the low-energy theorem <|14[1 dictates the zero- momentum limit of the physical correlator to be of the order of the 
light current quark masses, i.e. about an order of magnitude smaller. The additional screening contribution from Eq. 
(|172|l cancels most of the direct instanton contributions, however, and thereby restores consistency. The underlying 
mechanism becomes most transparent in the chiral limit where the low-energy theorem requires the subtraction term 
to vanish and the cancellation to become complete (cf. Eq. (|14|l ). Indeed, the zero-quark-mass limit of the screening 
contribution, 



n <T> (0) = = (16?r ^° )2 = 2Vn (174) 



(obtained from the definitions of F v and F v i and the mass relation Ijl68|l with m v — niK = <p>- = 0), exactly cancels the 
direct instanton contribution (|173|) and thereby restores consistency with the low-energy theorem. The topological 
charge is totally screened in the chiral limit since the massless Goldstone boson exchange generates infinite-range 
interactions. The presence of the screening contributions thus explains how the direct-instanton induced subtraction 
terms can be of equal size in both spin-0 channels whereas the low-energy theorems (|10|l and l|14|) require the size of 
their phenomenological counterparts to differ by an order of magnitude. This provides compelling evidence for the 
vital role of topological charge screening in the IOPE of the ++ glueball correlator. 

Furthermore, the cancellation between the direct-instanton- and screening-induced subtraction terms suggests a 
strategy for implementing the IOPE scale p into the screening contributions. As shown above, restricting the size of 
the direct instantons to p < \i~ x implies the replacement of n in Eq. I|173|) by n^ir = with ( < 1 (cf. Eq. (|148|) ~). 
Compliance with the low-energy theorem then requires the same replacement in the screening contributions l|174|) 

and therefore in ITp ar ' (Q 2 )- This is physically reasonable since only instantons with p < /i -1 ~ O (A scr ) can take 
part in the short-range screening mechanism. A conceptually similar restriction has already been implied by using 
the pointlike (i.e. p = 0) approximation in deriving Eq. (|169fl . 

We conclude this section with a first look at the quantitative impact of the screening contributions on the h Borel 
moments. To this end, we contrast Fig. 2 (a) of the k = moment, based on the perturbative and instanton-induced 
Wilson coefficients as given in Eq. I|149l) . with the same plot but including the screening contributions, in Fig. 2 (b). 
The previous deficiencies, i.e. both the negative sign and the positive slope of the Borel moment in the sq > 7 — 8 
GeV 2 region, are clearly resolved. In particular, lZp t o now decays monotonically and approximately exponentially 
with t and thus contains a clear signal for a pseudoscalar glueball resonance (and the rj , see below) which will be 
analyzed quantitatively in Sec. IVII (We recall that the negative regions of the Borel moment at smaller so and r are 
an artefact of the spike distribution and will disappear when realistic instanton size distributions are implemented, 
cf. Fig. 4 below.) 

In summary, we have found strong evidence for the direct instanton contributions alone to be a seriously incomplete 
description of the hard, nonperturbative physics in the ^ glueball correlator. Sum rules based on this selective choice 
of contributions |42j . i.e. without the screening contributions, are in several ways inconsistent. Reliable sum rules can 
be obtained only if the complementary physics generated by topological charge screening is properly included. 



VI. QUANTITATIVE ANALYSIS OF BOREL MOMENTS AND SUM RULES 

We now assemble the various IOPE contributions - (|64|) . (|129|l and l|172|) obtained above into our theoretical 
prediction for the continuum-subtracted Borel moments (|27|l . 

TZ G , k (r; So ) = ft<$ (r, s Q ) + Tlffi (r; s ) + fc.P^f (r) , (175) 

which form the left-hand side of the Borel sum rules (|28|l . These moments summarize the microscopic (i.e. quark- 
gluon level) information contained in the IOPE and the local duality continuum. Some additional input of a more 
implicit nature is needed, however, before quantitative glueball properties can be extracted from them by means of a 
sum rule analysis: a specific parametrization of the resonance sector (i.e. the number of isolated resonances and their 
shapes and widths), the determination of the duality threshold, the establishment of the "fiducial" r-domain in which 
all underlying approximations are expected to be reliable, and finally the choice of sum-rule optimization criteria. 

Before tackling these issues in Sec. IVIDI we will make an effort to obtain as much qualitative insight as possible 
from the IOPE moments alone. The main focus will be on pertinent features of the r- and so _ dependence as well as 
the scales which govern it, and an analysis of the various subtraction constants and their impact. The gained insights 
will help both to select the optimal strategy for the sum-rule analysis and to better understand the physics which 
underlies its predictions. 
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A. Input parameters and scales 

Before embarking on the quantitative analysis, we have to fix the values of the various constants and scales which 
appear in the IOPE and its moments i|175|) . All of them will be standard (or in their standard range). The dominant 
soft scale is set by the gluon condensate, 

(a s G 2 ) = (a s G a ^G atlu ) ~ 0.055 GeV 4 (176) 

(this value lies in the middle of the phenomenologically acceptable range (a s G 2 ) ~ 0.035-0.075 GeV 4 0123), whi ch 
determines the higher-dimensional gluon condensates as 

(gG 3 ) = (gfabcG^GfG^) ~ -1.5 (a s G 2 ) 3 ' 2 , (177) 

(a 2 s G*) s ~ ^(a s G 2 ) 2 , <a 2 G 4 ) p ~ f (a s G 2 f . (178) 

All condensate values refer to the renormalization scale fi ~ 0.5 GeV. The lattice estimate for the three-gluon 
condensate [z^ is sometimes replaced by the single-instanton estimate (gG 3 ) ~ 0.27 GeV 2 (a s G 2 ) . (Adopting different 
values of the lowest-dimensional condensates (inside their standard range) would affect mainly the predictions of the 
pseudoscalar k = —1,0 sum rules where the relative impact of the power corrections is largest.) Furthermore, the 
phenomenological subtraction constant in the pseudoscalar channel, determined by the low-energy theorem 114(1 . 
contains the quark condensate 

(qq) ~ - (0.24 GeV) 3 . (179) 

For the remaining IOPE parameters we adopt the values 

m u ~m d ~ 0.005 GeV (180) 
A Q c7Z5 = 0.2GeV (181) 
/i = 0.5 GeV (182) 

For the leading moments of the instanton size distribution we use the canonical values |l0f 

n ~ - fm~ 4 = 7. 53 x 10~ 4 GeV 4 , (183) 

p ~ - frn = 1.69 GeV -1 (184) 
3 

which also completely determine the finite-width size distributions discussed in Section IIV Al Our quantitative 
predictions will be based on the probably most realistic of these parametrizations, the Gaussian-tail distribution 185f) 
with the soft cutoff l|143fl and (3 — 3 GeV. The corresponding direct-instanton fraction £, defined in Eq. i|148|) . is 

- i />oo 

^ ^ = _ / {p) ^ a _ j ^ Q _ 66 (185) 

n n J 

which implies that about two thirds of all instantons are direct, i.e. small enough to affect the Wilson coefficients. (The 
(■-values from the exponential-tail distribution (|76|l and/or the alternative cutoff function ()144fl differ by about 1%.) 
Note that both the screening- and direct-instanton-induced subtraction constants are proportional to the instanton 
density n and therefore equally attenuated by renormalization (recall that F^^i, m 2 , <x fi), 

r4 scr) (o) -> o4r r) (o) , (186) 

n^W^cn^Co). (187) 

Finally, we recall the masses and couplings 

F 2 , = 0.543 GeV 6 , F 2 = 0.0886 GeV 6 , (188) 
m v = 0.55 GeV, m v > = 0.96 GeV (189) 

which determine the topological charge screening contributions l|172|) . 
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B. Qualitative behavior of the Borel moments 

In Figs. 3 and 4 we plot all four IOPE Borel moments (|175fl in both spin-0 glueball channels. The direct instanton 
contributions were calculated on the basis of the Gaussian-tail size distribution (|85|l and renormalized according to 
Eq. I|142fl with (3 = 3. These Borel moments will be the theoretical input for our sum rule analysis in Sec. IVI El 
For instructive purposes, i.e. to exhibit several qualitative features more prominently, the plots cover a region of the 
t — So plane which extends beyond the limits in which the sum-rule analysis is reliable. 

An obvious feature of all Borel moments is that they are monotonically increasing with so and monotonically 
decreasing with r in their physically sensible domains. This behavior stabilizes multi-parameter sum rule fits since 
the minimization routine is unlikely to get trapped in a local minimum of the deviation measure (cf. Sec. IVI Dfl . 
Furthermore, the monotonic increase of the r-slope with so, especially at small r, ensures that the largest resonance 
mass squared remains below the duality continuum threshold, as it should be. 

Another general property of all Borel moments is their increasing So-independence for r 3> Sq . This is a conse- 
quence of the Laplace suppression factor exp (— rs) in the dispersion integrals (|27() (which therefore have practically no 
support at s ^ so) and the damping of the instanton contributions at large s by realistic instanton size distributions 
(cf. Sec. HVEl . 

All Borel moments show an approximately exponential decay with increasing r (for so > 6 GeV 2 in the lowest 
h moments, which includes the physical region). Thus, all of them provide resonance signals in the sum rules l|28|l . 
After normalizing the moments to a common scale, the T-slopes of the h Borel moments are consistently larger 
than those of their ++ counterparts in the region r < 0.8 GeV -2 where the heaviest isolated resonance dominates. 
Thus t he glue ball mass scale is considerably larger in the h than in the ++ channel, in agreement with lattice 
results ,28, 43]. The same conclusion can be drawn independently from the ratios T^cfc+i/T^G.fc of adjacent moments 
in each channel. For k > and if one pole dominates, they are about equal to the square of the resonance mass and 
indeed consistently larger in the h channel. 

The pseudoscalar moments (especially those with k > 0) flatten out rather suddenly in the r-direction for r > 1 
GeV -2 . The visibly slower decay in the large-r region is a qualitative indication for the appearance of a second res- 
onance with a clear mass scale separation. The quantitative sum rule analysis below will confirm that the screening 
contributions have produced an additional signal for the 7]', as expected from "unquenching" the instanton contribu- 
tions. 

Several qualitative features of the Borel moments can be traced to the impact of finite-width distributions and 
renormalization on the direct instanton contributions. To exhibit them, we compare the k = moments of Figs. 
3 and 4 with those obtained from the traditional spike distribution, plotted in Fig. 5. Clearly, the neglect of 
interference between contributions from instantons of similar size in the spike distribution distorts the moments at 
small r: positivity and monotonic rise with so disappear. Moreover, the negative r-slope of 7Zs,o at small r turns 
positive at larger so and creates a "mountain ridge" . In addition, IZs.o itself turns negative in the latter region 
while IZpfi becomes negative mostly in the small-So (and small-r) region. In both channels, the maximal r-slope and 
consequently an upper bound on the glueball mass predictions is reached at intermediate so- 

The positivity violations generated by the spike distribution can be traced to the large-s oscillations in the instanton- 
induced imaginary part. As discussed in Sec. IIV El these oscillations affect mostly the small-r behavior of the moments 
and are a conceptually worrisome artefact. Their practical impact is moderate, however, since the positivity violations 
occur mostly outside of the fiducial r- and so-domains in which the sum rule analysis takes place. This is in marked 
contrast to the impact of the positivity violations which we have discussed in Sec. |Vl i.e. those which arise in the 
pseudoscalar Borel moments when the screening contributions are ignored. Indeed, a glance at Fig. 2 (a) shows that 
these occur mostly at intermediate and large r, i.e. exactly in the region where the sum rules are matched. They 
would make a meaningful sum rule analysis impossible and once more underline the necessity of the topological charge 
screening contributions. 

The renormalization of the direct instanton contributions generally increases the size of the moments at small r 
(due to the enhanced imaginary part at large s, cf. Fig. 1 (b)) and reduces it at the upper end of the fiducial 
r-domain (due to the effectively reduced instanton density n — > £n) ■ Both effects tend to increase the resonance mass 
predictions of the IOPE sum rules (cf. Eq. (|28|l '). Another effect of the improved treatment of the nonperturbative 
Wilson coefficients is a substantial reduction in the overall size of the direct instanton contributions. In addition, the 
perturbative contributions are enhanced by the 3-loop corrections (and additionally by the now favored, larger values 
for Aqcd, cf. Sec. Iffl A|l . Hence the previous dominance of the instanton- induced coefficients, found in Ref. 01 on 
the basis of the two-loop radiative corrections and the spike distribution, is considerably weakened. As a side effect, 
the derivation of the scaling relations between glueball and instanton properties [13J is obscured since the latter relied 
on the dominance of the instanton contributions at intermediate and large r. Nevertheless, these relations remain 
suggestive, especially because they are consistent with large- A^ c counting. 

A semi-quantitative upper bound on the sum-rule predictions for the glueball masses can be obtained from the 
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moment ratios considered above. Its derivation starts from the expressions 



T^G,k (r; s ) 
^ G , fc _i(r;s )+4 > on^ ) (0) 



^ (r;s ) = J- rV"?'<^ ( 19 °) 



(k = — 2) which are obtained by setting fc2 = in the sum rules l|28[) and solving for mci- (An alternative approach 
would be to take logarithmic r-derivatives of the Borel moments.) Upper bounds for the glueball masses are then 

found by searching for the position in the (r, So) plane where uiq is least sensitive to variations in both r and sq. 
These are either extrema or inflection points of Eq. (|190fl . Estimates of this type are often used instead of a full 
sum-rule matching analysis although they can accommodate only one resonance pole and are of limited reliability. To 

(2) 

give a specific example of this approach, we will use m s to establish a bound on the scalar glueball mass. From our 

Borel moments we obtain mi 2 ' as plotted in Fig. 6. One reads off an extremum (maximum) in t at t* ~ 0.2 GeV~ 2 
for So < 6.5 GeV 2 and an inflection point in so at Sq ~ 4.5 GeV. (In the analysis based on purely perturbative Wilson 
coefficients, one finds instead an inflection point of rrig in r and an extremum (minimum) in so Q.) Together, they 
yield the upper bound 

m s < mf (r*; s* ) ~ 1.7 GeV (191) 

for the scalar glueball mass. The analogous estimates from m^' 1 ^ result in somewhat smaller bounds. Nevertheless, 
these bounds should not be considered as a substitute for the sum rule results. Indeed, our quantitative sum-rule 
analysis in Sec. IVI El will show that they overestimate the ++ glueball mass prediction by about 35%. 

(2) 

Figure 6 furthermore reveals that m s has a stronger r-dependence than the analogous expression for the spike 
size distribution (cf. Ref. 13], Fig. 3). This might indicate that a one-resonance sum-rule analysis is somewhat less 
favored if realistic instanton size distributions and renormalization of the instanton-induced coefficients are taken into 
account. 

Finally, the Figs. 3 and 4 show that the k = — 1 Borel moments have a r-independent offset which becomes visible 
at large r. It is rather large and negative in TZs,-i while smaller and positive in TZp,-±. These offsets are due to the 
subtraction terms which the nonperturbative IOPE coefficients generate. (For this reason, Hq-i < does not imply 
positivity violations.) Their match to the subtraction constants on the phenomenological side of the sum rules is an 
important consistency criterion which we are going to discuss in the next section. 



C. Subtraction constants 



By design, the k = — 1 Borel moment included first-principle information provided by the low-energy theorem 1)10(1 
into the ++ glueball sum-rule analysis As pointed out in Sec. Ill Al an equally useful low-energy theorem l|14|) 
exists in the pseudoscalar channel and suggests to analyze the analogous h sum rule, based on the moment 7£p._i, 
as well. In order to prepare for this analysis, the present section investigates the conceptual and quantitative impact 
of the involved subtraction constants. 

Since the perturbative UV contributions to the subtraction constants are removed by definition in the low-energy 
theorems (cf. Sec. Ill AJ) and by renormalization in the perturbative IOPE coefficients, it remains to clarify the role 
and treatment of the nonperturbative contributions. The direct-instanton induced subtraction terms 

n s/p (°) = ±2Vfw ~ ±C x 0.95 GeV 4 (192) 

are part of the k = — 1 Borel moments (|105|l and (|129|) . Their size is significant but smaller than the size of the 
r-dependent contributions, in particular those due to direct instantons (at small and intermediate t). Moreover, for 
C = 1 (i.e. for the spike distribution) the instanton-induced subtraction constant in the scalar channel is more the 
50% larger than (and of the same sign as) the phenomenological LET value (cf. Eq. (|1U|) 

n (ph) (Q) = 327T ^ o ^ Gev4 ^ (ig3) 

Oo 

It would remain larger even when the largest available value for the gluon condensate, (a s G 2 ) ~ 0.07 GeV 4 , is used. In 
the pseudoscalar channel the discrepancy is yet more pronounced: the value (|192[) is more than an order of magnitude 
larger than the LET value (T2J| 

U { f ] (0) = (8tt) 2 mumd (qq) ~ -0.022 GeV 4 (194) 
m u + nid 
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(and again of the same sign). 

The impact on the k = —1 sum rules results from the fact that TZq _i has to be fitted to the decaying resonance 
exponentials and the "phenomenological" subtraction constants l|193(l and (jl94[l . As a general rule, the smaller the 
difference between the IOPE subtraction constants and the LET values, which makes up the remaining offset, the 
better will be the fit quality to the resonances (which is the only intrinsic reliability measure for QCD sum rules). 
The impact of the remaining imbalance can be rather subtle since it is largest towards the upper boundary of the 
fiducial r-domain where the match to the exponential resonance contributions becomes delicate. 

A glance at the above scales confirms that such a match is ruled out in the ^ channel, with its extreme discrepancy 
between LET and direct-instanton induced subtraction constants, if the nonperturbative IOPE coefficients arise 
exclusively from direct instantons. As we have shown in Sec. IV CI however, this inconsistency is overcome by the 
crucial topological charge screening correction 

rTp Scr) (0) = C(^f + ^f|~Cx 0.89 GeV 4 , (195) 

which cancels most of the direct-instanton contribution l|192|) and brings the total IOPE subtraction constant in line 
with the small LET value i|194[) . Previous analyses of h sum rules have discarded the k = — 1 moment and therefore 
missed valuable first-principle information from the low-energy theorem as well as a useful consistency check. 

One might wonder what happens to the rather delicate balance between the subtraction constants in the 
Nf — limit, i.e. in pure gauge theory or in the quenched approximation. In this case the above cancellation 
does not work since topological charge screening disappears. As is well-known, however, the LET (|14|> and thus the 
phenomenological value for the zero-momentum correlator is strongly affected by the absence of light quarks, too. 
Indeed, it becomes 

jjCph.gn) (Q) = _ (g7r) 2 ^(gn) „ _ Q m (1Q6) 

where we have used the standard (quenched) lattice value Y + qn \ — (0.18 GeV) 4 for the topological susceptibility [67I . 
which is in agreement with the Witten-Veneziano formula |8l|. Thus we find again perfect consistency: the LET 
subtraction constant becomes much larger and cancels the direct-instanton contributions by itself (assuming that the 
latter are not strongly affected by quenching). 

A special situation arises if the difference between the IOPE and LET values of the subtraction constant becomes 
so strongly negative that even the r-dependent contributions cannot prevent 

n G ^ 1 (r; SQ )+U^ h) (0) (197) 

from turning negative inside the fiducial r-domain |l!2fl . In this case the sum-rule fit to the (always positive) resonance 
exponentials in Eq. (|28[l is substantially worsened even if the negative offset in (I197|l remains relatively small. Exactly 
this situation is encountered in the scalar k — — 1 sum rule since 

Tl ( P h) (0) - n|j /+/) (0) ~ (0.61 - 0.95 C) GeV 4 < 0. (198) 

For £ = 1 (spike distribution) the imbalance between the r-independent terms is of about the same size with and 
without the direct instanton contribution. The main effect of the instanton contribution is to turn the sign negative, 
which causes the mentioned decline in sum-rule consistency. One option to deal with this problem is to assume 
that unreliably calculated soft contributions dominate the instanton-induced subtraction constant (in particular when 
using the spike distribution), as discussed in Sec. IIV CI and therefore to discard it completely. This strategy was 
adopted in Ref. [13j where the constant Ijl92|l was removed from the Borel moments (|105(l and (I129|l . It strongly 
improves the sum-rule quality and maintains the crucial mutual consistency with the predictions of the k > sum 
rules and with the LET (jl0|l (which is badly violated if nonperturbative contributions to the Wilson coefficients are 
ignored) This shows that the main reason for achieving a large glueball mass scale and LET consistency in 

the 1Zs,-i sum rule is the large overall size and slope of the r-dependent direct instanton contributions (and not the 
relatively small, instanton-induced subtraction constant, as suggested in Ref. [8fj|). 

Of course, discarding the instanton-induced subtraction constant completely is a very crude way of "renormaliza- 
tion". Our implementation of the IOPE renormalization scale on the basis of realistic size distributions accounts more 
accurately for the soft physics to be removed. Comparison of the above scales shows that this significantly improves 
the consistency with the LET: the difference between instanton- and LET-induced subtraction terms is reduced by 
C < 1, and for realistic values £ ~ 2/3 (cf. Eq. I|185|) ') the sum l|197|) remains positive over the whole fiducial r-domain 
Note that this renormalization maintains consistency with the LET in the pseudoscalar channel as well, while 
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simply discarding the instanton-induced subtraction constant would leave the compensating screening corrections out 
of balance (cf. Sec. IV (J() . 

In the context of this section it might be useful to recall that the large "phenomenological" subtraction constant had 
caused serious problems in the ++ glueball sum rule as long as the nonperturbative Wilson coefficients were ignored: 
the ensuing, weaker decay of the sum (|197l) generated a much smaller ++ glueball mass "prediction" (well below 1 
GeV) which was inconsistent with the predictions well beyond 1 GeV from the k > sum rules. It was pointed out 
in Ref. ^| that the missing, strongly decaying r-dependence can only reside in the Wilson coefficients and that it 
should be of nonperturbative origin, suggesting (together with other indications) direct instantons as its main source. 
The inclusion of the direct instanton contributions indeed overcame the mutual inconsistencies among the ++ sum 
rule predictions and restored the consistency with the LET. 

In summary, the direct-instanton induced subtraction terms and their renormalization play a rather complex role in 
the k — —1 glueball sum rules. In both channels they are essential for achieving consistency with the underlying low- 
energy theorems. In the pseudoscalar sum rules, this additionally requires strong cancellations with the indispensable 
topological charge screening contributions. These cancellations explain, in particular, how the equal size of the 
instanton-induced subtraction constants in both spin-0 glueball channels can be reconciled with the conspicuous 
difference between the LET values. 



D. Sum rule analysis setup 

In order to extract glueball properties from the IOPE Borel moments by means of a QCD sum-rule analysis, the 
parametrization of the phenomenological spectral functions and the matching criteria have to be specified. Of course, 
the complexity of the phenomenological side (measured by the number of parameters to be predicted) is limited by the 
information content and resolution power of the truncated short-distance expansion. Therefore, judicious decisions 
are required about the amount of detail to include, e.g. about the number of isolated resonances, their individual 
widths and shapes, specific assumptions on quarkonium admixtures or even explicit multi-hadron continua (beyond 
the local-duality approximation). Guided by the asymptotic nature of the spacelike IOPE (i.e. the factorial and 
not Borel-summable growth of the higher Wilson coefficients |8^ |) and its tr unca tion, we restrict ourselves in the 
following analysis to at most two resonance poles in zero-width approximation |l!4j and the local-duality continuum, 
as anticipated in Eq. (|28() . 

Any quantitative sum-rule analysis requires a numerical measure 6 for the deviation between both sides over the 
discretized fiducial r-domain (whose boundaries T min ,r max will be determined below). The iterative minimization of 
S up to the desired accuracy is then performed numerically. We will adopt the Belyaev-Ioffe measure |33( 




( 



fin (ft<g e > {n) - 8 k ^ h) {0),HG, k (rO) + £ , 



(199) 



with N = 100 grid points and r, = r min + i (r max — r min ) /N |115| . The constant £ is an offset to be added if otherwise 
the argument of the logarithm would become negative. An important requirement on reliable sum rule fits is that they 
should be stable, i.e. that the resulting hadron properties should be - inside the other typical errors of the analysis - 
independent of the starting values for the iterative minimization of 6. As mentioned above, the monotonic behavior 
of the IOPE Borel moments generally improves stability by reducing the likelihood for less than optimal local minima 
of 5. We have tested several alternative expressions for 6 (with, e.g., different weights for the deviations in the large- 
and small-r regions) and found them to change predictions of stable fits maximally at the one-percent level. 

The "fiducial" r-domain, in which the sum-rule analysis takes place, is designed to optimally exploit the physical 
information in the IOPE without leaving the region of validity of the involved approximations. Hence, one seeks the 
maximal r-interval in which the sum rules can be expected to be both reliable and predictive. Towards small r, the 
duality continuum in the Borel moments i|25|) increasingly dominates the glueball signal. In order to ensure that the 
sum rules remain sensitive to the glueball properties, we therefore fix r m j n (k, sq) by the standard requirement that 
the continuum contributions to the Borel moments must not exceed the resonance contributions, i.e. 



sp{cont) t _ \ 
T^G,k (T m i n ; So) 



0.5. (200) 



Generally, the higher Borel moments require larger values of r m i n . The above criterion therefore also helps to assure 
consistency among sum rule predictions from different moments. Furthermore, it typically reduces the fiducial domain 
of the higher-fc sum rules and thereby renders their fits somewhat less stable. (The results for sq, in particular, become 
less reliable with increasing k since the so-dependence of the Borel moments is largest at small r.) 
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The standard sum-rule criterion for determining the upper limit of the fiducial domain, r max , is to restrict the 
contribution of the highest-dimensional operator to maximally 10% of the total OPE contribution, i.e. 

C { 8 G) (r max ; M )(6 8 ) 

Ed=o,4,6,8 C d (r max ; n) l^Od J 

The neglected contributions from d > 8 operators should therefore remain small up to the onset of the asymptotic 
IOPE region. The above criterion is more stringent at smaller k where the condensates contributions have a relatively 
larger impact. (The size of the topological charge screening contributions also decreases with increasing k.) 

In the presence of nonperturbative IOPE coefficients, the criterion for r max requires some additional thought. Even 
for exclusively perturbative Wilson coefficients, the condition 12011) it is not too restrictive in the spin-0 glueball 
channels since the power corrections are unusually small. When the large direct instanton contributions to the 
unit operator are added, it can become almost ineffective. (In the pseudoscalar channel their impact is reduced by 
cancellations with the screening contributions and the larger so-values.) Another constraint on r max , arising from 
the requirement that multi-direct-instanton corrections to the Wilson coefficient should remain negligible, will then 
become more stringent. In order to formulate this criterion quantitatively, we adopt the rough estimate 

r < r max ^\{R~ 2p) 2 ~ 1.5 GeV- 2 . (202) 

For the higher Borel moments this criterion is typically more restrictive than (|201|l since the relative impact of the 
power corrections at large r decreases with increasing k. In view of its approximate nature, however, it is reassuring 
that the sum rule results are rather insensitive to variations in r max . (The sensitivity to T m - ln is larger, as expected.) 

Under the phenomenological parameters to be determined by the sum rule analysis, the duality threshold sq plays 
a special role. It is not associated with a glueball property but rather corresponds roughly to the squared mass gap 
between ground state and first radially excited state. Indications for a delayed onset of duality [8j| suggest, however, 
that this rule of thumb is invalid in the scalar glueball channel. (These indications could be tested on the lattice when 
reliable unquenched glueball spectra become available.) Hence, our only roboust expectation for s is that it should 
be larger than the squared mass of the highest-lying isolated resonance. 

Below, wc will determine s together with the glueball parameters from the sum rule fits, which turns out to be 
possible even in the presence of two isolated resonances. Nevertheless, it might be instructive to briefly comment on 
alternative strategies for obtaining the duality threshold. One alternative would be to divide the fitting procedure 
into two steps, by first constraining y/so to exceed the largest resonance mass by a constant amount As, 

y/s^ = m G + As (203) 

and by subsequently minimizing 6 as a function of As. Such a constant splitting is probably not an unreasonable 
assumption since the r-slopes of the IOPE Borel moments at small r increases rather monotonically with i/sq (cf . Figs 
3 and 4) in the sum-rule relevant r-region. (The couplings fs are much less so-dependent.) Although this two-step 
procedure can accelerate the minimization procedure (in particular in the two-resonance case), we did not find it 
necessary even in the pseudoscalar channel. Another strategy, adopted e.g. in Ref. |36j |. determines so by analyzing 
a related finite-energy sum rule (FESR) which roughly expresses a duality constraint. In view of the likely delayed 
onset of local duality this procedure might be misleading, however, in the scalar glueball channel. In Refs. fA Ht|. 
furthermore, so is required to render specific combinations of Borel moments minimally r-sensitive. In Ref. Q, finally, 
an upper limit on the FESR value for so is obtained by locating the extrema or inflection points of moment ratios 

like (jnoj. 

E. Results and discussion 

After having done the groundwork, we now proceed to the quantitative sum-rule analysis. This amounts to matching 
the continuum-subtracted IOPE Borel moments (|175fl to either one or two isolated resonances and, if k = —1, to the 
phenomenological subtraction constant (cf. Eq. J2EJ0- The decision about how many resonances to include will be 
made in the scalar and pseudoscalar channels individually, based on a comparative quantitative analysis. Of course, 
the more flexible two-resonance parametrization is almost bound to reduce the nominal fit errors at least slightly. 
Since this comes at the price of two more parameters to be determined, however, we will resort to the two-resonance 
parametrization only if it leads to a clear improvement of the fits. This is a necessary requirement for stable and 
physically meaningful predictions of additional resonance parameters. 
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1. Scalar glueball 

The analysis of the ++ glueball Borel sum rules on the basis of the spike distribution (with perturbative coefficients 
up to O (a s )) showed that a one-pole fit could almost perfectly match the IOPE moments [LI- This result left little 
room for improvement and therefore no indication for the presence of a second low-lying resonance with strong coupling 
to the scalar gluonic interpolator J3J). 

An extensive numerical survey of all four improved ++ Borel sum rules, based on the IOPE 11751) . reveals that 
they can be well fitted by only one ++ glueball pole, too. Again, there is no conclusive evidence for the presence 
of another low-lying, strongly coupled ++ resonance. A sufficiently exhaustive analysis of the scalar sum rules can 
therefore be based on the traditional one-pole plus duality continuum parametrization. The main part of the following 
discussion will deal with the results of this "benchmark" analysis. Nevertheless, there seems to be some indication for 
additional low-lying strength, as foreshadowed in the analysis of the moment ratio in Sec. EIH We will come back 
to this issue at the end of this section. 
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r min (GeV- 2 ) 


r max (GeV- 2 ) 
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1.06 
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0.8 


1.5 


1.22 


1.04 


1.63 


1.63 


2 


1.0 


1.5 


1.39 


1.01 


1.82 


2.13 



TABLE I: The fiducial domain, fitting error and predictions of the scalar glueball sum rule based on the k-th Borel moment. 

The numerical results of the one-resonance analysis of all four Borel-moment sum rules are collected in Tab. [I] 
They can be summarized in the overall prediction 

ms = 1.25 ± 0.2 GeV, fs = 1.05 ± 0.1 GeV (204) 

for the scalar glueball mass and coupling. The error assignment includes the estimated uncertainties in the input 
parameter values. The corresponding sum rule fits, in their individual fiducial r-domains, are shown in Fig. 7. 

Clearly, the overall fit quality is more than satisfactory. As anticipated, there seems to be no need for a second 
isolated (and rather narrow) resonance. The 3-parameter fits are inside typical errors independent of the starting 
values and also quite independent of the upper border of the r-domain. (Nevertheless, far-off starting values for 
the threshold sq can significantly slow down the minimization procedure since the So-dependence is rather shallow.) 
The instanton continuum contributions decisively improve the individual and mutual consistency of the sum rules 
(including the one associated with 1Zs,~i) and of their results. As expected, we find little sensitivity to details of the 
instanton size distribution as long as its overall scales are kept fixed. The sum rule consistency noticeably worsens, 
however, when finite-width distributions are used without subtracting the large-p contributions. 

The Hs,-i sum rule, singled out by the presence of subtraction terms (cf. Sec. IVI CI) . has played a notorious 
role in previous glueball sum rule analyses and deserves some more specific comments. Since the phenomenological 
subtraction constant is of significant size, the contributions of the perturbative Wilson coefficients alone were unable 
to generate a glueball mass of more than a few hundred MeV. As explained before, the direct instanton contributions 
with their large slope (cf. Fig. 7) resolve this problem 13]. The instanton-induced subtraction term, however, 
is overestimated by the spike approximation to the instanton size distribution: it so strongly overcompensates the 
phenomenological subtr actio n constant that the sum rule becomes inconsistent (the matching error increases by about 
2 orders of magnitude) |ll6fl . As anticipated in Sec. IVI CI this problem is resolved by realistic p-distributions and 
renormalization of the instanton contributions. Indeed, Fig. 7 shows that the fit quality of the improved 1Zs,-i sum 
rule matches that of its higher- moment counterparts, and all sum rules yield mutually consistent predictions for the 
glueball properties. 

The impact of the subtraction constants on the sum rule analysis increases towards larger r since the dominant 
r-dependent contributions decrease. The sum-rule fit in Fig. 7 shows that the residual discrepancy between LET- and 
IOPE-induced subtraction terms can be comfortably compensated by the small condensate contributions which are 
enhanced in the k — — 1 sum rule. Even larger discrepancies, encountered e.g. when using smaller values of the gluon 
condensate, could be accommodated. The k — 1 sum rule would then predict a larger value of sq and a marginally 
larger glueball mass. 

In addition to the sum-rule fits for the total IOPE Borel moments, Fig. 7 also shows the contributions from 
perturbative and nonperturbative Wilson coefficients separately. Comparison with the analogous figures in Ref. [l3| 
(based on the spike distribution) confirms the discussion of Sec. IVI Bl the perturbative contributions are somewhat 
enhanced by the O (a^) corrections while the direct instanton contributions are significantly reduced in the improved 
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IOPE. (The duality continuum threshold decreases somewhat, too.) As a result, perturbative and nonperturbative 
coefficients are now of comparable size. Moreover, the slopes of the instanton contributions are reduced. Both effects 
manifests themselves in a 20% smaller prediction for the scalar glueball mass. However, the about threefold increase 
of the glueball decay constant fs predicted in Ref. 01 (relative to previous sum rule analyses) remains intact. This 
implies that the reduced size of the Borel moments on the left-hand side of the sum rules l|28ll is compensated by the 
increased perturbative contributions and, in particular, by the smaller glueball-mass factors on the right-hand side. 

Our large result for fs predicts strongly increased partial widths for the radiative decay of the heavy quarkonia 
J/ip and T (and others) into scalar glueballs and therefore has relevance for experimental glueball searches |85j. In 
particular, it allows for an improved analysis of the existing T — > 7/0 decay data of the CLEO collaboration |86| 
and the forthcoming larger samples from CLEO-III. It will also be interesting to compare our prediction to the first 
calculation of fs on the (quenched) lattice which is in progress [87j ■ Since fs is related to the glueball wave function 
(or Bethe-Salpeter amplitude) at the origin, our large value predicts a strongly concentrated wave function and thus 
an unusually small size of the scalar glueball. Similar results have been found in the instanton liquid model |15J and 
on the lattice [l4j . 

Our prediction for the central value of the scalar glueball mass is about 10-20% smaller than the results of our 
previous analysis based on the spike distribution, and lies 10-40% below the quenched lattice results m^' = 1.4—1.8 
GeV |2S|. (All raw lattice data agree within statistical errors (about 40 MeV). The much larger range quoted above 
reflects the ambiguities in setting the mass scale in the absence of direct experimental input (cf. the talk of Bali under 
Rcf. 43]).) A similar reduction (25%) of the ++ glueball mass was found in unquenched lattice simulations (on 
still rather small lattices with very limited statistics and quark masses of about 70 MeV) although this might 
be dominantly a lattice artifact. On general grounds, however, one would expect light-quark effects and quarkonium 
admixtures to lower the quenched masses at least somewhat. Due to its smaller size, furthermore, the scalar glueball 
should be particularly susceptible to the momentum dependence of the sea-quark induced vacuum polarization which 
modifies the color-dielectric properties of the vacuum at short distances |88| . Since this part of the vacuum polarization 
is neglected in the quenched approximation, one would expect the quenched mass predictions to be less reliable in 
the scalar channel. This expectation is supported, e.g., by the recent analysis of the radial-excitation and Regge 
trajectories of the known isoscalar mesons |8J|, including those newly established by the Crystal Barrel Collaboration 
in proton-antiproton annihilation. The standard ++ glueball candidates above 1 GeV, i.e. the fo (1300), fo (1500) 
and fo (1750) resonances, are found to lie solidly on qq trajectories and to fit well into a flavor nonet classification. 
Instead, the iVmatrix analysis predicts a light and broad scalar glueball state in the 1200-1600 MeV region |89| . 
compatible with our result. An even lighter (and similarly broad) glueball state, centered around 1 GeV or a bit 
larger, is expected in mixing schemes which assume only one ++ multiplet below 1.8 GeV |90j . 

The fc-dependence of our predictions for the scalar glueball properties in Tab. [i] shows a certain systematics. While 
the results for the coupling are practically fc- independent (within sum-rule accuracy), the predictions for the mass 
increase with fc for fc > 0. Since fs is associated with the integrated strength of the glueball, it enters all sum rules 
(|28|l in the same (i.e. fc-independent) power. The fc-dependence of fs in Tab. U thus gives an idea of the typical 
uncertainties of the matching analysis. The variations in the mass predictions, on the other hand, are larger and 
systematically increase with fc. This suggest that the glueball strength is distributed over a rather broad s-region: 
since the higher moments weight the spectral function more strongly at larger s, they will then predict a larger 
pole mass (cf. Eq. (|28|) ~). (The k = — 1 sum rule result is not conclusive in this regard since it receives additional 
contributions from the subtraction terms.) Our predicted mass range can therefore be regarded as a rough lower 
bound on the width of the scalar glueball, i.e. Ts > 0.3 GeV, similar to the width found in the iVmatrix analysis 

m. 

The gap As between glueball mass and continuum threshold, as defined in Eq. I]203p . is relatively constant among 
the fc > sum-rule results, 

As = 0.5 ± 0.1 GeV, (205) 

while it is about twice as large in the fc = — 1 sum rule (As = 0.94 GeV), due to the subtraction term. The assumption 
of a common shift As for all sum rules in a simplified analysis (cf. Sec. IjVI Dfl l would therefore fail in the lowest 
moment. The relatively early onset of the continuum may indicate, incidentally, that the first excited scalar glueball 
state lies around 2 GeV. Although the quenched lattice spectrum predicts this excitation well beyond 2 GeV and 
beyond the lowest ^ glueball state, its mass might again be lowered by quarkonium admixtures. 

We had argued above that the quality of the one-pole sum rule fits leaves little room for physically significant 
improvement and provides no clear evidence for an additional low-lying resonance with strong coupling to the gluonic 
interpolating field. Nevertheless, a residual r-dependence of the Borel moment ratio m s was found in Sec. IVIBI and 
the fit quality of the Hs.o sum rule is somewhat lower than that of the other ones (cf. Tab. [I]and Fig. 7). Since this 
sum rule is probably most sensitive to low-lying strength, in particular at large r where it would be hidden in 1Zs,-i 
by the subtraction constants, it is tempting to attribute its reduced fit quality to some broad, low-lying structure 
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missing on the phenomenological side. Amusingly, this would be consistent with the iVmatrix analysis of the ++ 
meson data j8^ which tentatively identifies the broad "<r resonance" /o (600) as a probably gluon-rich exotic state. 

Unfortunately a two-resonance analysis of the lZs,o sum rule does not help to settle this issue: due to the high 
quality of the one-pole fits, it can hardly reduce the deviation measure further and thus becomes unstable. (If 
nevertheless performed, it seems to favor an only somewhat smaller mass msi ~ 0.8 — 1.2 GeV of the lowest-lying 
resonance and a heavier glueball at ms2 ~ 1.8 GeV.) We conclude that significant low-lying ++ strength with a 
sufficiently strong coupling to the gluonic interpolator @, if it exists, is probably to broadly distributed to be resolved 
by a sum rules analysis. 

2. Pseudoscalar glueball 

As in the scalar channel, we start by surveying the quantitative behavior of all four sum rules in order to determine 
how many isolated resonances are required. The result is opposite to that in the scalar channel: all sum rules, 
especially those derived from the lower Borel moments, clearly favor two resonances, with a large separation between 
their masses. The improvement over the one-pole fits is substantial. In fact, the latter do not only produce significantly 
larger errors but also tend to destabilize the sum rules since either the stronger decay at small r or the weaker one at 
large r, but not both, can be matched to one resonance. 

The emergence of a second, relatively low-lying isolated resonance is of course expected. In fact, its traces were 
already visible in the qualitative decay behavior of the Borel moments (cf. Sec. IVI B|l . Moreover, our discussion in 
Sec. anticipated the emergence of large rj' intermediate-state contributions to the h glueball correlator, as a 
consequence of the anomaly-induced ^'-coupling to the topological charge. The results of the quantitative sum rule 
analysis, listed in Table |nj indeed confirm this expectation: the central mass value of the lighter resonance reproduces 
the rj' mass. Our central value for the coupling f n ' is somewhat larger than its phenomenological value 0.82 GeV 78], 
likely because of Ty-admixtures. (Note that we treat f v ' on the same footing as the glueball decay constant, i.e. we do 
not extract the conventional factor 1/ */2iV/.) In any case, the result for f v ' is probably our least accurate sum rule 
prediction since it originates from the small residue of the large-r tail. 

A two-resonance analysis, with five independent hadronic parameters to be determined, generally stretches the sum- 
rule resolution to its limits. In fact, there is no guarantee that such an analysis will be stable. The large separation 
between the resonance masses in the channel, however, improves the situation decisively since it assigns mutually 
almost exclusive roles to the two poles: the heavier h glueball has to fit the small-r region and consequently 
decays so fast that it cannot significantly "contaminate" the large-r tail, which is mostly generated by the rj' . This 
scenario is corroborated by the fact that the predicted rj properties are, in contrast to the glueball properties, almost 
so-independent. A glance at Fig. 4 shows that the small-r behavior of the moments indeed varies much more strongly 
with so than the tails. It turns out that the clear mass separation renders the five-parameter fits stable and makes a 
quantitative sum rule analysis possible. 

As another consequence of the large gap between the resonance masses, the relative strength of the rj and glueball 
signals, oc (my /mp) 4+2fe , decreases strongly with k since the Borel moments weight the large-s region of the spectral 
function by a factor s k . For this reason, the glueball predictions of the higher-moment sum rules become more stable 
and the quality of one-pole fits improves with k, reflecting the diminishing impact of the rj . Nevertheless, the one-pole 
analysis remains inherently unstable since the resonance exponential can either match the IOPE moments at large or 
small r, but not over the whole fiducial region. Hence the two-pole fits continue to be superior even for k — 2. The 
impact of the condensate contributions, incidentally, also decreases with k since the higher ^-moments are obtained 
by derivatives with respect to — r and thus decrease the size of the power corrections at large r. 



k 


r min (GeV- 2 ) 


r max (GeV- 2 ) 


m„/ (GeV) 


U (GeV) 


m P (GeV) 


fp (GeV) 


so (GeV 2 ) 


5 x 10- 3 


-1 


0.2 


1.2 


0.97 


0.88 


2.12 


0.43 


6.79 


3.22 





0.25 


1.3 


0.81 


1.31 


2.32 


0.76 


9.14 


3.89 


1 


0.55 


1.4 


1.05 


0.84 


2.08 


0.72 


6.63 


1.59 


2 


0.55 


1.5 


1.08 


1.10 


2.20 


0.79 


7.31 


4.86 



TABLE II: The fiducial domain, fitting error and predictions of the pseudoscalar glueball sum rule based on the k-th Borel 
moment. 

The four h sum rule fits are displayed in Fig. 8 and produce the results contained in Tab. [HI Their central 



39 



values yield the predictions 

rap = 2.2 ±0.2 GeV, f P = 0.6 ± 0.25 GeV, (206) 
rry = 0.95 ± 0.15 GeV, = 1.05 ± 0.25 GeV, (207) 

again including an additional error due to estimated input parameter uncertainties. (One should keep in mind that 
the rj 1 pole probably receives some strength from ^-admixtures.) As in the scalar channel, the figures also show the 
contributions from the perturbative and nonperturbative IOPE coefficients separately. 

In the graph for lZp t -i, we have additionally plotted the topological charge screening contribution alone (dotted 
line). Its negative sign is a consequence of the screening-induced subtraction term. Nevertheless, the total nonper- 
turbative unit-operator coefficient (dashed line) is positive, due to the (positive) instanton-induced subtraction term. 
The balance between these subtraction constants is instrumental in reconciling the k = — 1 sum rule with the low- 
energy theorem (JT2J, as argued in Sec. IVI Gl Without the topological charge screening contributions this sum rule 
would obviously be inconsistent and impossible to fit. The dependence on the LET parameters (quark masses and 
condensate), incidentally, is much weaker than in the scalar k = — 1 sum rule since the phenomenological subtraction 
constant is an order of magnitude smaller. The cancellation between the individually much larger IOPE contributions 
is therefore a highly nontrivial consistency requirement of the anomalous axial Ward identity (cf. Sec. Iv|) . 

In the higher moment sum rules, the screening contributions remain essential and generate, besides the rf resonance, 
a clear h glueball signal. The plots of the nonperturbative Wilson coefficients in Fig. 8 contain the cancellations 
with the direct instanton contributions and demonstrate why it would be detrimental to ignore the topological charge 
screening contributions (cf. Sec. IV! The screening contributions affect the behavior of the moments mostly 
at t > 0.5 where their impact becomes comparable to that of the perturbative and direct instanton contributions. 
Positivity violations due to the direct instantons and the dissolution of the glueball signal would manifest themselves 
in this r- region if the screening contributions were ignored |42| . 

Despite the larger values of so, the size of the direct- instanton contributions is not much smaller than in the scalar 
channel. Our discussion of the finite-width distribution and renormalization effects in the previous section applies 
to the most part here, too. However, due to the cancellations among the nonperturbative contributions their overall 
impact on the pseudoscalar sum rules is much more moderate. (The conventional criterion (|201|l therefore determines 
Tmax only in the k = 2 sum rule.) As a consequence, one expects the results to become closer to those of the older h 
glueball sum rule analyses which neglected nonperturbative Wilson coefficients altogether. Indeed, Ref. 36] found 
772 p = 2.3 ± 0.2 GeV which is compatible with our result. 

Our prediction for the glueball pole residue, however, is about twice as large as the value f P — 0.30 ± 0.05 GeV 
of Ref. |36| . This enhancement is due to the remaining nonperturbative contributions, the higher-order perturbative 
corrections and the additional rj' pole. Since fp is related to the scale of the light-cone distribution amplitude for the 
pseudoscalar glueball, our result can be used, e.g., to determine the 77 — > Gpir cross section at large momentum 
transfer [9l| ■ As in the scalar channel, the increased prediction for the coupling f P and the consequently larger partial 
width of radiative J/ip decays into pseudoscalar glueballs js^j are of relevance for experimental glueball searches. 

Our prediction for the h glueball mass lies inside the range of quenched lattice results, m,p = 2.1 — 2.5 GeV 
|28j . (The range of values again reflects scale-setting ambiguities.) Unquenched simulations are still at an exploratory 
stage, with correspondingly large errors for the higher- lying glueball masses |43j. The pseudoscalar mass, however, 
is close to that for the lowest tensor glueball and in the same range as the quenched results. The fact that our 
prediction for the h glueball mass is close to the quenched lattice results while that for the ++ glueball mass is 
significantly smaller may be related to the smaller ++ glueball size. Since sea-quark-induced vacuum polarization, 
which is missing in the quenched approximation, affects mainly the short-distance properties of the glueball wave 
functionals, one might expect the quenched mass predictions to be subject to larger dynamical-quark corrections in 
the scalar channel. 

The if-matrix analysis of Ref. |89| cannot identify narrow-resonance candidates for a ^ glueball in the above 
mass range since both 77 (1295) and 77 (1440) are found to lie on linear quark-antiquark trajectories. The Gaussian 
sum rule analysis of gives an about 20% larger 0~ + mass than ours, outside of the range of lattice results. Due to 
the absence of the crucial screening contributions the underlying IOPE is inconsistent with the axial Ward identity, 
however, and these results cannot be trusted. 

In contrast to the results in the scalar channel, the fc-dependence of the pseudoscalar glueball mass shows no 
particular systematics, although the coupling is again practically /c-independent for k > 0. It might be tempting 
to speculate that this implies a less homogeneous distribution of the h glueball strength, perhaps generated by 
two (or more) neighboring resonances. In any case, it seems that no useful estimate for (or bound on) the width of 
the pseudoscalar glueball can be delineated from the variation of the pseudoscalar glueball properties with k. (The 
monotonic increase in the predictions for the 77' mass with k, as well as the larger fluctuations in the coupling, are 
probably due to contaminations by the increasingly dominant glueball signal.) 
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The derivative of the topological charge correlator at Q 2 = 0, plays an important role in the analysis of the proton 
spin content and the related structure function |34[ . The strong cancellations among the nonperturbative contributions 
to the pseudoscalar IOPE might explain why QCD sum rule estimates of x' on the basis of purely perturbative Wilson 
coefficients seem to be sufficiently stable ;93j. Nevertheless, this analysis should be repeated on the basis of the full 
IOPE. 



VII. SUMMARY AND CONCLUSIONS 

The central themes of this paper are the derivation of nonperturbative short-distance contributions to the spin-0 
glueball correlators and a comprehensive sum-rule analysis of their predictions for glueball properties. The dominant 
nonperturbative contributions to the operator product expansion of these correlators turn out to be (semi-) hard, 
i.e. to reside in the Wilson coefficients. Their soft counterparts, contained in the condensates, play a comparatively 
minor role. Both sources of hard nonperturbative physics considered in this paper, direct instantons and topological 
charge screening, are associated with the topology of the vacuum gluon fields or, equivalently, of the QCD gauge 
group. The main benefit of analyzing the hard nonperturbative contributions by means of a short-distance expansion 
is that it allows an analytical and largely model- independent treatment. The complementary bulk information on the 
soft physics, which enters through a few condensates and the instanton size distribution, can be straightforwardly 
imported from other sources. 

Direct instanton contributions to other hadron correlators were calculated previously, and several important effects 
with diverse physical impact were found in the corresponding QCD sum rules (including those for the ++ glueball). 
However, their evaluation was restrained by three major approximations: (i) all instantons were taken to be of the 
same size, (ii) the renormalization of the instanton-induced Wilson coefficients was ignored and (iii) the contributions 
of only the instanton nearest to the interpolator arguments were taken into account explicitly. The first of these 
approximations becomes exact in instanton vacuum models at large N c and the third can be regarded as the leading 
term in an expansion in the instanton density. The second, although convenient, is not strictly associated with a 
systematic expansion. 

An important point on our agenda was to assess the validity of these approximations and to improve upon them. We 
have removed restriction (i) by developing and implementing minimal expressions for realistic, finite-width instanton 
size distributions which incorporate the currently available information on its scales, overall shape and limiting 
behavior. These expressions are fully determined by the standard values of the instanton density and the mean 
instanton size. To improve upon approximation (ii), we have devised a simple renormalization procedure which 
explicitly restricts the contributing instantons to the direct ones, i.e. to those which are smaller than the inverse 
operator scale. These developments are channel independent and will be useful in other hadron correlators as well. 

The improvements in the direct-instanton sector have several notable effects on the glueball sum rule analysis and 
its results. Finite-width distributions considerably slow the asymptotic decay of the instanton contributions at large 
momentum transfer, which reduces the glueball mass predictions. They furthermore resolve artefacts in the IOPE 
which previously distorted the Borel moments and contaminated the results for the scalar glueball. As an additional 
benefit, realistic instanton size distributions tend to enlarge the fiducial r-domain of the sum rule analysis. The 
gauge-invariant renormalization of the instanton-induced Wilson coefficients significantly improves the consistency 
with the underlying low-energy theorem in the scalar channel, and sum-rule consistency in general. Renormalization 
also reduces the overall size of the instanton contributions. The reanalysis of our previous IOPE sum rules in the 
scalar glueball channel shows that the beneficial impact of the direct instanton contributions is enhanced by the above 
improvements. The consistency among different Borel moments and with the low-energy theorem is consolidated and 
the previously deficient and usually discarded lowest-moment sum rule becomes one of the most stable and reliable. 

In the pseudoscalar channel the emerging pattern is more complex: the sum rules with purely perturbative Wilson 
are found to be consistent with the corresponding low-energy theorem (which had not been appreciated before since 
the pertinent lZp_i sum rule was not analyzed). However, this consistency is lost entirely when direct instanton 
corrections are included and additional physics, both (semi-) hard and nonperturbative, is needed to restore it. 

Our third improvement of the IOPE consists in identifying this missing physics as due to topological charge screening 
and in implementing the screening contributions. This is our conceptually farthest-reaching extension of the OPE 
since it introduces for the first time nonperturbative physics beyond single instantons into the IOPE coefficients. This 
physics turns out to be highly channel-selective: it almost exclusively affects the h glueball correlator (and the 
related 7/ correlator) which is proportional to the topological charge correlator and thus specifically tuned to those 
light-quark induced correlations which produce topological charge screening. Under the suggestive assumption that 
the topological vacuum charge is mostly due to instantons, the implementation of its screening can be regarded as 
an improvement upon approximation (iii) mentioned above: the small instanton packing fraction does not suppress 
the exceptionally strong and short-ranged correlations between opposite topological charges as strongly as others. 
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Therefore, in this specific case the nearest-instanton approximation is insufficient and the screening corrections have 
to be added in the 0~ + glucball (and rf) channel. 

We have found compelling evidence for the screening contributions to be an indispensable complement to the direct 
instanton contributions. They restore consistency with the axial Ward identity and thereby overcome, in particular, 
the above-mentioned problem with the low-energy theorem. This manifests itself e.g. in the screening contributions 
to the subtraction constant in the lowest-moment sum rule which cancel most of the direct instanton contributions. In 
the chiral limit, this cancellation becomes exact. It also explains in a natural way how the equal size of the instanton- 
induced subtraction constants in both spin-0 glueball channels can be reconciled with the order-of-magnitude difference 
in the sizes of their phenomenological counterparts. (In the absence of light quarks this difference would practically 
disappear together with the screening contributions, so that the sum rules would remain intact.) 

The topological charge screening corrections lead to the emergence of a strong rf resonance signal and necessitate a 
two-resonance analysis of the h sum rules. Remarkably, the large gap between the rf and h glueball masses allows 
for a simultaneous prediction of all associated resonance and threshold parameters. The impact of the direct instanton 
contributions differs strongly in both spin-0 glueball channels. In addition to their opposite sign (a consequence of the 
instanton's self-dual field strength), the larger continuum threshold and the strong cancellations with the screening 
corrections drastically modify the role of the instantons in the h channel and counterbalance their repulsion. In fact, 
ignoring the screening contributions would, besides violating the anomalous Ward identity, lead to spectral positivity 
violations and the disappearance of the h glucball signal. Pscudoscalar glueball sum rules with unscreened direct 
instanton contributions are therefore invalid. 

Our comprehensive numerical analysis of all eight Borel sum rules in both spin-0 glueball channels reveals a rather 
diverse pattern of glueball properties. The hard topological physics in the IOPE coefficients turns out to strongly 
affect the sum rule results and to generate several new predictions. In the scalar channel, the improved treatment 
of the direct-instanton sector reduces our earlier (spike-distribution based) result for the ++ glueball mass by about 
20%, to tylq — 1.25 ± 0.2 GeV. Although still consistent within errors, our new central mass prediction is smaller 
than the quenched lattice result. However, light-quark effects and especially quarkonium admixtures are expected 
to reduce the quenched masses, and the first unquenched simulations indeed show a tendency towards smaller scalar 
glueball masses. Moreover, our mass prediction is consistent with the broad glueball state found in a recent if -matrix 
analysis of the scalar meson spectrum which includes the new states recently identified in the Crystal Barrel data. 
The systematics in our results from different Borel-moment sum rules indicates a rather large width of the scalar 
glueball, Ts > 0.3 GeV. A similarly large width is found in if-matrix and mixing analyses and could also be expected 
from unquenched lattice simulations with realistic quark masses, due to the increased number of open decay channels. 

Our prediction for the glueball decay constant, fs — 1.05 ± 0.1 GeV, is several times larger than the value 
obtained when ignoring the nonperturbative Wilson coefficients. This result implies an exceptionally small glueball 
size, in agreement with some lattice and instanton liquid model evidence. The strong concentration of the ++ 
glueball (Bethe-Salpeter) wave function is therefore at least partially explained by the strong instanton-induced 
attraction between gluons in the scalar channel. Moreover, its small size makes the scalar glueball more susceptible 
to the momentum-dependence of the color-dielectric constant arising from sea quarks. Since this part of the vacuum 
polarization is missing in the quenched approximation, one would expect larger errors in the quenched predictions for 
the ++ glucball properties than for their ^ counterparts. This might explain why our mass prediction deviates 
more strongly from the quenched results in the scalar channel. Furthermore, our enhanced prediction for fg implies 
substantially larger partial widths of radiative J/ip and T decays into scalar glueballs. It is therefore of importance 
for experimental glueball searches, in particular for the interpretation of the recent CLEO and forthcoming CLEO-III 
data on T — > 7/0 and other decay branches, and for related measurements of scalar glueball properties. 

Although hints of some gluon-rich spectral strength in the ++ channel well below 1 GeV can be detected in the 
sum rule analysis, we find no clear evidence for a low-lying and sufficiently narrow glueball resonance in the region of 
the /o (600). (The 1Zs,-i sum rule would "predict" a spurious low- lying resonance if the direct instanton contributions 
were ignored. However, this result is inconsistent with both the low-energy theorem and the higher-moment sum rules 
and therefore obsolete.) 

In the pscudoscalar glueball channel, the topological charge screening contributions do not only resurrect the sum 
rules but also have a strong impact on their quantitative predictions. Due to the cancellations with the direct 
instanton contributions, the overall size of the nonperturbative IOPE coefficients remains below about 20% of their 
perturbative counterparts and the detrimental problems encountered when ignoring the screening contributions are 
resolved. Although the origin of the individual nonperturbative contributions is more diverse and complex than in the 
scalar channel, their overall impact is therefore smaller. This is one of the reasons for the considerably weaker binding 
of the h glueball. Nevertheless, the hard nonperturbative contributions modify qualitative features of the 0~ + Borel 
moments to which the matching analysis is particularly sensitive, and they are vital for achieving consistency among all 
moment sum rules and with the axial anomaly. The strong cancellations among the nonperturbative contributions to 
the pseudoscalar IOPE may also explain why QCD sum rule estimates of the derivative of the topological susceptibility, 
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which plays an important role in the analysis of the proton spin content, seem to be sufficiently consistent without 
nonperturbative Wilson coefficients. Their reanalysis on the basis of the full IOPE is in progress. 

The quantitative sum rule analysis results in the values nip = 2.2 ± 1.5 GeV for the pseudoscalar glueball mass and 
fp = 0.6 ± 0.2 GeV for the decay constant. Our mass prediction lies inside the range obtained from quenched and 
unquenched lattice data. The coupling is again enhanced by the nonperturbative Wilson coefficients, but less strongly 
than in the scalar channel. The consequently larger partial width of radiative quarkonium decays into pseudoscalar 
glueballs and the enhanced 77 — > Gpir cross section at high momentum transfers will be relevant for the experimental 
identification of the lowest-lying h glueball and help in measuring its properties. 

Our extended IOPE should be useful for the calculation of other spin-0 glueball properties as well. Quantitative 
estimates of the already mentioned production rates in gluon-rich channels (including J/ip and T decays) and char- 
acteristic glueball decay properties and signatures, including 77 couplings, OZI suppression and branching fractions 
incompatible with qq decay, would be particularly interesting. 
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APPENDIX A: BOREL MOMENTS FROM PERTURB ATIVE WILSON COEFFICIENTS 

In this appendix we list the contributions from the perturbative OPE coefficients to the Borel moments Ck (t, sq) 
with k £ {—1, 0, 1, 2} for both spin-0 glueball correlators (before RG improvement). In terms of the continuum factors 
Pk and the exponential integral E\ , defined as 
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APPENDIX B: INSTANTON INTEGRALS 



An explicit expression for II 



c+') 



can be obtained by Fourier transforming its dispersive representation (|118f) : 
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as anticipated in Eq. (|92|l . The hypergeometric function 2-Fi (a, b, c, z) [54 , l57j is denned as the analytical continuation 
of Gaufi' hypergeometric series (except if c is a nonpositive integer — n and neither b nor c equal an integer —m with 
to < n) 
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which imply in particular, 
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The integral IjBll) can alternatively be done by several other methods, e.g., by introducing Feynman parameters, 
which leads to 101 
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but is generally less convenient for practical calculations and for the study of analyticity properties. 

A straightforward way to calculate the Fourier transform of Hq + ^ (x 2 ) starts from the original integral 
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one then immediately obtains 
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In order to perform the Borel transform of l|B13(l . it is convenient to start from an integral representation f |57|. 
Eq. 8.486.15) for the McDonald function, 
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The three remaining parameter integrals are elementary. Their evaluation leads to 
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The higher moments with k > — 1 follow, as previously, by differentiation with respect to — r, 
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Figure captions: 

Fig. 1: Dircct-instanton induced imaginary part of the ++ glueball correlator, ImIIg /+ ^ (— s), obtained on the basis 
of (a) the spike and (b) finite-width instanton size distributions. In (b) the results for the exponential-tail 
distribution without (dotted) and with (dashed) large-p cutoff as well as for the Gaussian-tail distribution with 
large-p cutoff (full line) are plotted. Note the difference in scale between (a) and (b). 

Fig. 2: The k = continuum-subtracted IOPE Borel moment of the 0~ + glueball correlator, calculated on the ba- 
sis of the spike distribution (a) without (left panel) and (b) with (right panel) topological charge screening 
contributions. 

Fig. 3: The continuum-subtracted IOPE Borel moments of the scalar glueball correlator (calculated on the basis of the 
Gaussian-tail distribution, and renormalized at the operator scale), as a function of t (x-axis) and So (y-axis). 
All units are appropriate powers of GeV. 

Fig. 4: The continuum-subtracted IOPE Borel moments of the pscudoscalar glueball correlator (calculated on the basis 
of the Gaussian-tail distribution, and renormalized at the operator scale), as a function of r (x-axis) and s 
(y-axis). All units are appropriate powers of GeV. 

Fig. 5: The k = continuum-subtracted IOPE Borel moments of the ++ (left panel) and h (right panel) glueball 
correlators, as obtained from the spike distribution. 

Fig. 6: The square root of the ratio between the k = 2 and k — 1 continuum-subtracted IOPE Borel moments of the 
scalar glueball correlator. 

Fig. 7: The individual contributions to the optimized Borel sum rules in the scalar glueball channel: the continuum- 
subtracted IOPE moments (full line), the ++ glueball pole (and subtraction constant, for k = —1) contributions 
(bullets), the direct instanton contributions (dashed) and the contributions of the perturbative Wilson coeffi- 
cients (dash-dotted). 

Fig. 8: The individual contributions to the optimized Borel sum rules in the pseudoscalar glueball channel: the 
continuum-subtracted IOPE moments (full line), the h glueball and tf (and subtraction constant, for k = — 1) 
contributions (bullets), the contributions of the nonperturbative Wilson coefficients (due to both direct instan- 
tons and topological charge screening) (dashed) and the contributions of the perturbative Wilson coefficients 
(dash-dotted). The figure for k — — 1 additionally shows the screening contributions by themselves (dotted). 
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gauge fields live. 

RG-invariance of Os to all orders in a s can be achieved with the help of the QCD /3 function by defining Os (x) = 
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[97] In Ref. |2C| a somewhat heuristic attempt was made to generalize the above type of low-energy theorem to the pseudoscalar 
sector in pure Yang- Mills theory. Since we are working within full QCD, whose three light quark flavors alter the behavior 
of the _+ correlator drastically, we have no use for these results here. 
[98] In studies of the OPE at large orders d, instantons of a different type were employed as a tool for assessing the impact of 
near-origin singularities on the asymptotic OPE behavior |35| . (In the context of an operational, OPE-based definition of 
quark-hadron duality such effects are referred to as "local duality violations".) 
[99] In Ref. |3q| it has been argued that the results of the standard-OPE sum rules (without direct instanton contributions) 
for the pseudoscalar glueball for Nf = is close to that of real QCD (Nf = 3) because the leading perturbative Wilson 
coefficient has a weak ^/-dependence (weaker than in the scalar channel) and since mixing with quarkonium seems to be 
smaller in the h than in the ++ sector |37l| . 

[100] By now there exists a rather large body of work on direct instanton contributions to IOPE and Borel sum rules in classical 
hadron channels [24J which uncovered important new instanton effects in hadron structure (in the nucleon channel e.g. the 
stabilization of the chirally-odd nucleon sum rule |lal . the emergence of a new, stable sum rule for the nucleon magnetic 
moments, the reconciliaton of sum-rule and chiral perturbation theory predictions for the neutron-proton mass difference, 
a nonperturbative dynamical mechanism of isospin violation |44T| etc.). 

[101] In the IOPE of the pseudoscalar meson correlator, analogous instanton-enhanced power corrections are associated with 
chiral-symmetry breaking operators and therefore turn out to be important I55j . 

[102] This conclusion is supported by quark and gluon condensate estimates on the basis of the one-instanton approximation 
which are not far from the phenomenological values, and by the fact that their multi-instanton corrections (evaluated in 
the ILM) are typically of the order of 10-20%. (I thank Sasha Dorokhov for interesting discussions on this point.) 

[103] if IR divergencies are encountered, see e.g. Refs. (6l| 

[104] Note that the large-p cutoff affects higher-O (h) corrections as well since it effectively restricts the internal momenta of 
the quantum fluctuations. 

[105] Another way to see that the large-p cutoff cannot entirely exclude soft contributions derives from the singularity structure 
of the direct-instanton contributions 19211 at x = 0, as discussed in Section II V Bl Indeed, potential instanton-induced 
power corrections are not removed by the large-p cutoff because the small- x singularity structure is not affected by 
large instantons (for physically sensible n(p)). (The nucleon correlator IOPE [l^| provides a simple example: instanton- 
induced power corrections are present despite the large-p "cutoff" implicit in the spike distribution.) Note, incidentally, 
that the self-duality I15H of the instanton's field strength severely restricts potential instanton-induced power corrections 
in quark-based correlators (to at most a few terms) |63|. 

[106] Additional /^-dependence of the condensates due to nonperturbative physics (e.g. due to instantons 64, 65]) can of course 
not be excluded. Our results and general sum-rule experience suggest, however, that it should be similarly weak for 
H ~ 2 - 3Aq C d ■ 

[107] As discussed in Section HVI multi-instanton effects are generally negligible, due to the small instanton packing fraction. 
In the pseudoscalar meson channel, where direct instanton contributions have exceptional strength |5^|. this was shown 
explicitly |46j . 

[108] As a consequence of the chiral anomaly (see below), the same holds for the correlator of the divergence of the flavor- 
singlet axial current (with the quantum numbers of the rjo) in the chiral limit. This correlator is thus equally sensitive to 
topological charge correlations. 

[109] Note the direct analogy with long-range Coulomb interactions which get screened into short-range Yukawa interactions 
between the electric charges of an electrolyte or a QED plasma, thereby turning the photon effectively massive. 

[110] At first, there seems to be a certain analogy with the effective quark mass which instanton-induced quark zero modes 
acquire due to interactions with ambient, soft vacuum fields. This effective mass also goes beyond the isolated-instanton 
approximation. However, it originates from soft subgraph corrections (of mean field type) to hard subgraphs and therefore 
does not affect the ^-dependence of the Wilson coefficients. The topological charge screening corrections, in contrast, have 
a very pronounced x- and channel dependence. 

[Ill] The calculation of the quark-mass dependence of gluon condensates on the lattice is impeded by small signal-to-noise ratios 
and problems with reaching physical light-quark masses. The quark-mass dependence is expected to be non-negligible. 

[112] When r — > oo this will eventually happen for any negative difference (and does, of course, not imply that the spectral 
function violates the positivity bound). 

[113] Although their numerical values are quite close, one should keep in mind that the physics entering Ilg + ^ (0) and FI^' 1 ' (0) 

is mostly complementary: as part of a Wilson coefficient, 11^ + ' (0) receives dominantly (semi-) hard contributions while 

Y[(p h ) (o) ig renormalized by subtracting the hard (perturbative) fluctuations and therefore dominated by soft modes (note 

the appearance of the gluon condensate in the LET IIOH 'I. 
[114] Finite resonance widths require exponential resolution and are therefore in principle inaccessible to the standard OPE |83|. 

The situation is more complex when including nonperturbative Wilson coefficients which themselves introduce exponential 

contributions. In any case, the zero-width approximation in QCD sum rules does not require the corresponding, physical 

resonances to be particularly narrow. 
[115] Due to the generally monotonic decrease of the glueball Borel moments with r, scale- invariant measures favor an improved 

matching in the large-r region. In two-resonance fits this favors the smaller-mass resonance. 
[116] The same should probably hold for the corresponding Gaussian sum rule which, perhaps for this reason, was excluded in 

the analysis of Ref. (l6|. 
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